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[1] DERIVE fiir den Mathematikunterricht, W. Kopf 

[2] Historische Aspekte im Mathematikunterricht, M. Kronfellner 

[3] Einfiihrung in den TI-89, Teil 1 und Teil 2, B. Kutzler, bk teachware 

[4] Computergestiitzer Physikunterricht, Experimente zur Mechanik - Auswertung mit dem CAS 
des TI-92, K-H. Keunecke, Texas Instruments 

[5] Introduction a la geometric avec la TI-92, J.-J. Dahan, ellipses, ISBN 2-7298-9877-8 

[6] Faire des Mathematiques au lycee avec des calculatrices symboliques, Luc Trouche, 

IREM des Montpellier, 34095 Montpellier Cedex 05, email: irem@math.univ-montp2.fr 


Tj\X Praxis, Raymond Seroul und Silvio Levy, Birkhauser, ISBN 3-7643-2823-1 

TgX - das geniale Satzsystem Donald Knuths - setzt sich auch bei nicht ganz professionellen Schrei- 

bem und Vielschreibem durch. Mit der aktuellen Version von Scientific Workplace gibt es nun auch 
eine Windows- Implementierung von TpX in Kombination MAPLE. 

Nach der englischen Ubersetzung liegt nun dieses franzosische Standardwerk fiber das Arbeiten mit 
TEX auch in einer sorgfaltig an den deutschen Sprachgebrauch angepassten Fassung vor. Das Buch 
richtet sich sowohl an den Einsteiger, wie auch an den fortgeschrittenen Anwender. 

Die Kapitel konnen selektiv gelesen werden, viele Beispiele sind in den Text eingestreut und verffih- 
ren den Lesenden zum Probieren. Insgesamt 16 Kapitel verteilen sich auf mehr als 400 Seiten. Die 
Kapitelfiberschriften reichen vom Zeichensatz von TpX fiber Seitenlayout zu Tabellen und Tabellie- 

ren. Dem wichtigen Formelsetzen sind fast 40 Seiten gewidmet. Ein ausffihrliches Glossar und Wor- 
terbuch runden das empfehlenswerte Buch ab. 


Symbolic Rewriting Techniques, Manual Bronstein a.o., Ed., Birkhauser, ISBN 3-7643-5901-3 

This volume 15 of the series "Progress in Computer Science and Applied Logic” contains invited and 
contributed papers to the Symbolic Rewriting Techniques workshop, which was held in Ascona, 1995. 
Symbolic rewriting techniques are methods for deriving consequences from systems of equations, and 
are of great use when investigating the structure of the solutions. 

Most of the contributions deal with Grobner bases e.g. The Computation of Grobner Bases Using an 
Alternative Algorithm (JApeT), Grobner Fans and Projective Schemes (D.MalT), etc. In total 14 con- 
tributions can be found within 288 pages, each of them with an extended reference list. I must confess 
that it is not an easy to read book and it seems not to be an introduction into this field of mathematical 
research, but it might be of high value for people who have some experience in this subject matter. 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English arti- 
cles will be warmly welcomed nonetheless. 
Your contributions will be edited but not 
assessed. By submitting articles the author 
gives his consent for reprinting it in DNL. 
The more contributions you will send, the 
more lively and richer in contents the DE- 
RIVE Newsletter will be. 


Next issue: March 1999 

Deadline: 15 February 1999 


Preview: Contributions for the next issues 

3D-Geometry, Reichel, AUT 

Graphic Integration, Linear Programming, Various Projections a. o., Bohm, AUT 
A Utility file for Complex Dynamic Systems, Lechner, AUT 
Examples for Statistics, Roeloffs, NL 

INTEGQ - a Package for the 77s, Propper & Bohm, GER & AUT 
Fractals and other Graphics, Koth, AUT 
Implicit Multivalue Bivariate Function 3D Plots, Biryukov, RUS 
Quaternion Algebra, Sirota, RUS 

Various Training Programs for Students on the TI-92, Bohm, AUT 
A Critical Comment on the "Delayed Assignation" := =, Kiimmel, GER 
Sand Dunes, River Meander and Elastica, The lighter side ...., Halprin, AUS 

Type checking, finite continued fractions Welke, GER 

On the Resolution of the Linear Differential Equation , Candel, ESP 

Share Holders’ Considerations using a CAS, Bohm, AUT 
Kaprekar’s “Self Numbers”, Finite Groups, Schorn, GER 
Graphic solution of Linear Programming Problems in the Tl, Kirmse, GER 
reversions of Spigot Algorithms for n and e, Witthinrich, GER 
Problems from Physics, Magiera, POL 

and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA , Koller, AUT; Meagher, AUT; 
Wiesenbauer, AUT; Aue, GER; Speck, NZL, and and and 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 

Herausgeber: Mag.Josef Bohm, Herstellung: Selbstverlag 


The DERIVE-NEWSLETTER is the Bulletin 
of the DERIVE User Group. It is published 
at least four times a year with a contents of 
44 pages minimum. The goals of the DNL 
are to enable the exchange of experiences 
made with DERIVE as well as to create a 
group to discuss the possibilities of new 
methodical and didactical manners in teach- 
ing mathematics. 

We have established a section dealing with 
the use of the 77-92 and we try to combine 
these modern technologies. 

Editor: Mag. Josef Bdhm 
A-3042 Wurmla 
D'Lust 1 
Austria 

Phone: 43-(0)660 3136 365 
email: nojo.boehm@pgv.at 
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Volker Loose, Germany 

Hi everyone, is it possible to have axes with logarithmic scale in DERIVED 
Peter van den Sanden , Netherlands 

There are no logarithmic scales in DERIVE. However, you can create logarithmic scales (using any 
base a) with a few substitutions which can be declared by functions: 

fOO := 


flog(x,a):=f(a) the function you want to plot 

f_log_y(x, a) := L0G(f0O, a) => logarithmic x-scale 

f_log_y(x, a) := L0G(f0O, a) =^> logarithmic y-scale 

x 

f_log_xy(x, a) := LOG(f(a ) , a) =^> logarithmic x- and y-scale 


0.3 

f 00 := 5-x 
f_log(x, 10) 

f_log_y(x, 10) 
f_log_xy(x, 10) 

1/5 

f(10000) = 50-10 

APPRQX(LGG(f (10000), 10)) = 1.S9S970004 



0.5 1 1.5 


2 2.5 


3 


3.5 4 4.5 
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Michael Johann. Germany 

I am using DfW 4.08. How can I make it to get only the real solution for x 3 = -1? I tried it by defining 
the range of x and Algebra Status - in vain. 

DNL\ There were several answers in the DERIVE - news pointing at the right choice of the 
branch for evaluating (-1 ) 1/3 . A nice answer was given by Terence Etchells, Liverpool: 

Hi Michael, 

Hmm, I’m intrigued! Why do you only want the real solution(s)? If you have a vector of solutions 
(e.g. [x = 1, x = -1, x = 1 + i, ...] ) and you wish to pick out the real ones, then a function 
REAL SOLVE does the trick: 


This request is obsolete for the recent version of DERIVE. Terence’s trick is no more neces- 
sary and it doesn’t work because the output of SOLVE has changed. Using DERIVE 6 it is 
easy to control the domain of the solutions. 


3 

SOLVE Cx = -1) = O = -1) 


3 

SOLVE Cx = -1, x) 


r 

x 


1 

2 


V X 

2 


1 4^ ■ i 

+ v x = -1 

2 2 


6 

SOLVE (x - 1 = 0) = (x = -1 v x = 1) 


6 

SOLVE O -1=0, x) 

1 1 73- i. 

X = - - V X = - + V X 

2 2 2 2 


3 

SOLUTIONSCx = -1, x, Real) = [-1] 


3 

SOLUTIONSCx = -1, x) = 


1 fi.L 

- 1 , + 

2 2 


6 

SOLUTIONSCx = 1, x, Real) = [1, -1] 

6 [ 1 Jl-L 

SOLUTIONSCx = 1, x) = 1,-1, + 

2 2 


1 fi.L 1 fi-L 

— - V X = + V X - -1 V X - 1 

2 2 2 2 


1 73.1.' 


2 2 


1 fi- t 1 Jl- L 1 73-d. " 

2 2 ' 2 2 ' 2 2 . 


DNL: Another FAQ (frequently asked question) appeared in the DERIVE-News followed by a 
very extended answer given by Johann Wiesenbauer, which seems to be interesting enough 
to be presented for you all: 

I was wondering if anyone knows a command for DERIVE that is similar to EVSUB from MuMath or 
SUBS in MAPLE. I want to evaluate f(x) = x 3 + 2x + 1 when x = 5 without using Manage Substitute. 
Can anyone help? 

Blain Priddle 

Some could. See Johann’s answer on the next page. Although we do have now the SUBST- 
command which is very powerful I don’t hesitate to reprint Johann’s proposals, Josef. 
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Now we can use the SUBST-command not only for substituting one variable. We can substi- 
tute a list of variables - and we can substitute for subexpressions: 

3 

SUBST (x + 2.x + 1, x, 5) = 136 
3 3 

SUBSTO -y + 2-x-y + 1, [x, y], [5, 3]) = 3406 


SUBST 


x + 


3 2 ^ 

, x , a 


2 

x + 2 


a + 2 


+ a 


But unfortunately this does not hold for all cases. As you can see, I was able to substitute for 
x 2 but not for (x+1 ) or (x+1 ) 2 . 


SUBST 


O + D + 


O+l) +2 


x + 1, a 


SUBST 


O + D + 


2 > 

•, 0+1) , a 


0 + 1 ) + 2 


x + 2«x + 3 
3 


+ a 


1 im 

X + 1 -+'EL 


1 "im 

x+a - 1 


O + D + 


O + D + 


■\ f 

= 1 im 

x + l-*a 


(x+1) + 2 J 


(x+1) + 2 j 


x + 2-x + 3 
3 

2 

x + 2«x + 3 


+ a 


2 ^, 


+ O + 1) 


3 2 

+ a 


2 

a + 2 


If you can do without the headline you can take the TABLE-command instead of Nurit’s func- 
tion from 1999. 



" 4 

S 

SO 

SOO 

TABLE([2.x, 20.x, 200.x], x, [4, 7, 25]) = 

7 

14 

140 

1400 


_ 25 

50 

500 

5000 


DNL: The next topic - Piecewise defined functions - appeared after a question posed by 
Allan S. Hugo. He wanted to define and then to plot a piecewise defined function like: 


y =/(*) = < 


x 0 < x < 1 

x 2 1 < x < 3 

2x + 3 x > 3 
undef. else 


I will summarize the very fundamental discussion which followed: 
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Andrei Jakobcic, Slovenia 

Try this: 

f (x) : =IF (x<0 , ?,IF(x<l, x, IF(x<3, x A 2, 2x + 3))) 


Johann Wiesenbauer, Vienna 

Of course, this is the appropriate way of solving this problem, though purist might prefer this version 
of f(x): 

f (x) : =IF (x>0 , IF(x<l, x, IF(x<3, x A 2, 2x + 3))) 


Ralph Freese, Honolulu 

You could also use DERIVED CHI-function. CHI (a, x, b) is 1 for x between a and b and 0 other- 
wise (assuming a < b). So graphing 

f (x) :=CHI (0,x, 1) *x + CHI (l,x, 3) *x A 2 + CHI ( 3 , x, inf ) * ( 2x+3 ) 

also does what you want. 

Johann Wiesenbauer, Vienna 

What looks very elegant at the first glimpse (most people don’t even know of the existence of the 
built-in CHI-function!) it is wrong at the second. The function 

f (x) : = CHI (0, x, 1 ) *x + CHI ( 1 , x, 3 ) *x A 2 + CHI ( 3 , x, inf ) * (2x+3 ) 

yields 0 for negative values of x, but should be undefined in this region. I would have liked this defini- 
tion, but I’m afraid, at least one IF is necessary e.g. 


f(x): = IF(x > 0, CHI(0,x,l)*x + CHI (l,x, 3) *x A 2 + CHI (3, x, inf ) * (2x+3) ) 


Due to the definition of CHI ( a , x , b ) for x = a and x = b, which is weird in my opinion, f(l), and f(3) 
yield a question mark instead of the correct values. (By the way, the same also goes for the version of 
Ralph Freese and my proposal for a correction.) Moreover, f(0) is undefined. Like it or lump it, the run 
off the mill solution using nested IF-statements (essentially due to Andrej Jakobic) seems to be the 
only correct one. 


f(>0 := x(0. 1)"X + x(l, X, x, 

-1 0 
0 0 
1 ? 


TABLE(f GO ■ x, -1, 5) = 


2 4 

3 ? 

4 11 

5 13 



(You can observe the gap for x = 3! Take care to 
edit the Chi-function in uppercase CHI, Josef.) 
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David Stenenga, Honolulu 

The function CHI (a, x, b) may not be that familiar to most DERIVE users but it is quite useful for 
certain types of problems. In our lab manual (see http://www.math.hawaii.edu/CalcLabBook/ [ * ] ), 
Ralph Freese and I have a nice example in connection with spline approximations. 


The idea is that you are given a(nx 2)-matrix which represents the xy-coordinates of some sampling 
of function values. We assume that the x- values are strictly increasing: x\ < x 2 < x 3 <... < x n but the 
y- values are arbitrary. The idea with spline functions is to define a piecewise polynomial function f(x), 
which is a polynomial f {x,k) on the interval x k < x < x k +i. Adjoining polynomials are subject to certain 
conditions such as continuity and perhaps some smoothness. 

Using nested IF-statements would be impractical unless the number of points is pretty small. Instead, 
we do the following (file PIECEWISE.MTH): 


# 1 : 

# 2 : 

#3: 

#4: 

# 5 : 

# 6 : 


[F(x , k) :=, data :=] 

DIMENSIGN(data) - 1 

I F(x, k)-x(data , x, data ) 

k=l k,l k + 1,1 


Example: Piecewise linear interpolation: 


Ffx, k) := 


data - data 

k + 1,2 k , 2 

data - data 

k + 1,1 k , 1 


> (x - data ) + data 

k , 1 k, 2 


Sample problem for piecewise linear "interpolation: 


fool 


data 


1 1 
2 0 


L 3 1 J 


#7 : Simplify #2 

(2 - x) ■ SIGN(x - 3) 

#S: + | x - 2 

2 


x — 1 1 + 


x 


2 



Of course, this technique will give a zero value for points outside the interval [xi, x n ] but then this can 
be corrected as noted by Johann. 

n The website is still active. You can find a file containing some DERIVE-tools at 
http://math.hawaii.edu/wordpress/category/undergrad-posts/computing-posts/derive-posts/ 
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Johann Wiesenbauer, Vienna 

Hi all, 

In general, the function CHI (a, x, b) is certainly very useful, but as we have seen its current defini- 
tion for x= a and x = b may cause problems, because its values are ambiguous, viz. 1/2 ± 1/2. Thus 
strictly speaking it’s not even a function! 

For this reason I’d like to suggest a more general DERIVE - function CHI (a, x, b, c, d) (hoping that 
Albert Rich is all ears) that is both a “genuine” function and also one that seems to live up to all (?) 
expectations: 

1. CHI (a, x, b, c, d) coincides with the current CHI ( a , x , b ) for a<x<b. 

Additionally, its values for x = a and x = b are c and d, respectively. In particular, this includes 
the indicator functions for the intervals [a, b\, (a,b\, [a,/?), (a,b) as well as the case c = d = ! / 2 , 
which plays an important role on some occasions, e.g. when dealing with Fourier series. 

2. CHI (a, x, b, c) isCHI(a,x,b,c,l-c) .This includes the indicator functions for the inter- 
vals (a,b] and [a, b) and again the case c = 1/2. These versions of CHI can be used to construct 
piecewise defined functions such as in David’s example, but the results will be “genuine” func- 
tions, i.e. they are uniquely defined for all x. When using nonlinear functions as in example 


f (x) : = CHI ( 0 , x, 1 ) *x + CHI ( 1 , x, 3 ) *x A 2 + CHI ( 3 , x, inf ) * (2x+3 ) 

this fundamental requirement is usually not fulfilled if x belongs to the set [xi, x 2 , x n ]. 
(It must be said though that for x\ or x n - or even both - an adjustment of the function value may 
be necessary. 

3. CHI (a, x, b) is defined according to its current definition in order to achieve total compatibil- 
ity. 

In my opinion, this modification of CHI is worth considering. What does the DERIVE community 
think of it? 


Cheers, Johann. 


It is interesting how the TIs are treating piecewise functions. In some respects they behave 
similar - and this is not really surprising but Nspire has a special feature. 


ifl^T FF? Y F]:t T FHt y ft Y Fe^ Y 'I 
| t i-— |fl 1 gebra Calc lother Pr gm 1 0 |c 1 ean Up | || 

1 1 I IP — n — TV 1 1 1 1 1 


IU2 

signCx) + 1 

'stepCx) 2 

1 step^x - a) - st-ep(x - b) + chi (a 7 x ? b^ 

Done 


... >*x A 2+chi<3 ..x,. »>^<2x+3>-^g<x>l 

Mil IN FiliD EKflCT FUNC H/50 




- #■:: Y FH 
Cel y |Header 

FmV-V 

1 

X 

yi 

y2 


-2.QQQQQQQ 

undef 

0.00000000 


-1.0000G0G 

undef 

0.00000000 


0.00000000 

0.00000000 

0.00000000 


1 . 00000000 

1 . 00000000 

signtG. 0G... 


2.G0G0G0GG 

4.G0G0G0GG 

4.G0G0G0GG 


3.0G0G0G0G 

9.0G0G0G0G 

9. 0G0G0G0... 


4.G0G0G0GG 

1 1 . 0G0G0GG 

1 1 . 0G0G0GG 


5.00000000 

13.0000000 

13.0000000 


lv2<x>=sign<G_ >+l 

m 

1 MAIN 

RAD EXACT 

FUNC 


DATA 






cl 

c2 

c3 | 

m 

c5 

c6 

1 

-1 

1 

G 

undef 



2 

0 

3 

G 

G 



3 

1 

5 

sign... 

1 



4 

2 

7 

4 

4 



5 

3 

9 

9*Cs... 

9 



6 

4 

11 

11 

11 



7 

5 

13 

13 

13 




c4=seg<f<cl[k]> J .k J .l J .10> 
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The Voyage 200 (and TI-92) does not know the CHI-statement. But following a contribution 
of our Canadian experts (DNL#91) it is easy to implement. Instead of IF in DERIVE we have 
to use the when-construct. When investigating the table we can observe the same deficien- 
cies of the CHI-defined function. Problems will appear when generating a value table in the 
Data-Matrix Editor. In column c3 I defined c3 = g(cl) - which works but f(cl) in column c2 
does not work correctly (f(1) = 1!!). So one cannot apply a when-generated function on a 
whole column. What to do? You can follow the construction of c4: define again a sequence 
addressing each single cell in column cl. 

Turning to Tl-Nspire we can do the same as with the Voyage 200 - and we face the same 
problems. 


/7(x):=when(x<0,undef,vvhen(x<l,x,when(x<3,x^,2- x+ 3 ))) 

Done 

( \ sign(x)+l 

stepyx ) : 

2 

Done 

chi(x ) : =step(x~b) -step {x~b) 

Done 

p(x):=x • chi(o,x, l)+x^* cfa’(l,x,3)+(2* x+3)* chi{3 7 x, 00 ) 

Done 



The left half of the screen shows the function value tables. Corresponding to the Data-Matrix 
Editor we want to generate a value table in the Lists & Spreadsheet Application. There we 
find that it is not possible to apply the “seq-trick” from above (see col. E next page). What we 
can do is presented in col. D: calculate the value in D1 and then copy down the cell in the 
column. This works!! 

But we have a better solution with Tl-Nspire: use the tablet for generating a piecewise de- 
fined function. This yields a perfect definition which meets all our requirements. 
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A 

B 

c 

D 

E 

- 


=f1(a[]) 

=f2(a[]) 


=seq(f1(a[k]), 

1 

-1 

1 

0 

#UNDEF 

piecewise(u... 

2 

0 

3 

0 

0 

piecewise(u... 

3 

1 

5 

±1+1 

1 

piecewise(u... 

4 

2 

7 

4 

4 

piecewise(u... 

5 

3 

9 

9*(±1+1) 

9 

piecewise(u... 

6 

4 

11 

11 

11 

piecewise(u... 

7 

5 

13 

13 

13 

oiecewiselu..^ 

<J^ )[> 

D1 

=/44 ] 

ii 


Use the tablet for defining the piecewise defined function: 


•§- □“ V5 75 e D hi> ft (ii <° {g 

|d| ™ [□□] [S] 1H „|p i D f D 

f? dn J adB ^ 



)• cfa'Opr, 00 ) 


/3(x): 


undef, x<0 
x, 0<x<l 

x 2 , l<x<3 

,2-x+3,x>3 


Done 

undef 

Done 
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Josef Bohm, Wtirmla 

Working with the TI-92 in my first form I made the following experience: I prepared an example for a 
written test: 


A linear system is given: 


ax + by = a + 2y-3x 
2a x - 3b y = b - 3y + 2x 


a) Solve this system for x and y. 

b) Find at least one pair (a,b) such that the system has no resolution. 

My idea was to check if the students were able to bring the system in a form to pick out the coeffi- 
cients for writing them into the respective matrix. Then they should apply rref() and interpret the 
solution in order to find values for a and b which let “disappear” the denominator. 

(As Tl-Nspire behaves pretty the same as TI-92 and Voyage 200 as well, I present the 
Nspire screens, Josef) 




One - female - student gave a very quick response. She wrote a = -3 and b = 2 without using the de- 
nominator and explained - during the test: Look at the first row: a = -3 and b = 2 leads to a contradic- 
tion, thus the system cannot have any solution!) I was happy to have such students! 


For the second group I just exchanged a and b on the left hand side in the system given above - with- 
out checking the further procedure at home. During the test I encountered after the good surprise de- 
scribed above a bad one, too: 
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The solution for x was a huge expression (the students could have taken the solution for y but they 
were too confused to do so). I encouraged the students to either do so or to copy the x-solution into 
their test forms - and made up my mind for the future to check even the “easiest” examples. 


Back home I applied factor on that crazy ex- 
pression for x and checked the same procedure 
with DERIVE using row_reduce(equb), wich 
worked properly - as expected. So I wrote a mail 
to David Stoutemyer reporting my “problem” and 
very soon I could receive his answer. 



David Stoutemyer, Honolulu 

This is a nice example of how reordering variables can make a big difference: Variable “ 0 ” is more 
main than variable “6”, and the choice about whether or not to combine rational expressions over a 
common denominator depends on the relative primeness of the two denominators. In this case the 
common denominator form is more compact, but the CAS doesn’t know that in advance and can’t 
afford to try all combinations. 

Aloha, David 

There was another question concerning systems of linear equations which initiated a nice 
discussion. Several answers followed. Til publish two of them which might be useful for you. 

Ingo Kronberger, 

I’d like to get a program for DERIVE with that I can do a Gaussian algorithm. I’d like the program to 
show all steps DERIVE is performing. Does someone have such a program or an idea how to do one? 
Thanks a lot. 

Johann Wiesenbauer, Vienna 

Here are some ideas that might be useful in this connection. Suppose you have a system of linear 
equations such as 

x + y + z = 6 
2x-y + z = 3 
3x + 2y-2z = l 

For my solution I need two functions: 

PIV0T(a, i, j, dir) := IF(dir, a, VECTOR(IF(di r .m_ < dir-i, a^m., IF(m_ = i, 
aj,i /aj,i j, j , aim_ - a 4 ,m_ 4 ,j/a 4 ,i jj - a^i)) , m_, DIMENSION(a)) , PIV0T(a, i, j, 1)) 

The pivot element a X] of the matrix A is used to force all elements of A of it (or above it, if di r = -1) to 
zero. Moreover the row a x is scaled in such a way that the new a X] becomes 1. If the fourth parameter 
di r is omitted, it is assumed to be 1 by default. 
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NEWA(u) := a := u 


This tiny auxiliary function is used for “self-assignments”, i.e. for assignments a := u where the ex- 
pression u contains a. The double assignment a:= = is supposed to do this, but it doesn’t yet work in 
DfW. 

Now a step-by-step solution could look like this: 


#4: 


#5: 


# 6 : 


#7: 


#S: 


1116 
2-113 
3 2-2 1 


NEWA(PIVOT (a, 1, 1)) = 


NEWA(PIVOT (a, 2, 2)) = 


1 1 

0 -3 
0 -1 

1 1 

0 1 
0 0 


1 6 
-1 -9 

-3 -17 

1 6 

1 

3 

3 

14 

-14 


NEWA(PIVOT (a, 3, 3, -1)) = 


NEWA(PIVOT (a, 2, 2, -1)) = 


110 3 
0 10 2 
0 0 13 
10 0 1 
0 10 2 
0 0 13 


This is just an outline of what could be done, but it should be easy to adapt the method above accord- 
ing to your own wishes, in particular, if you want smaller steps or other matrix operations such as 
swapping rows (cf. the utility-file VECTOR.MTH for this). 


NEWA is still working, the double assignment was not “reanimated” in later DERIVE- 
versions. There is one fact in Johann’s procedure which can be improved, it is necessary to 
take variable “a” for the matrix to be transformed. I changed this and made the input a little 
bit more comfortable, Josef. 


#9 : m 


nem(mx, i, j, dir) := 

#10: Prog 

m := PIVOT (mx , i, j , dir) 


r 4 1 1 6 20 


#12: newmCb, 1, 1) = m 



1 

4 

11 

2 


3 

3 

2 


0 -20 


#11: b := 


2-163 -10 

3 2 -2 S 0 


0 


3 

4 


11 7 

-13 

4 2 


L 4 -3 2 1 2 


0-6 1-3 -IS J 
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#13: newn(m, 2, 2) = m := 


#14: newnfm, 3, 3") = m := 


0 0 -21 -3 

1 1 


u u u 


#13: newmOn. 4. 4. -1) = m := 


3 


2 

3 

0 

40 

3 

21 

190 

11 

11 

3S6 

330S 

11 

11 

1 

o 

3446 

4 

193 

11 

o 

40 

3 

3 

1 0 

176 

193 

0 1 

1634 


0 10 0 


0 0 10 


0 0 0 1 - 


3490 

193 

192S 

193 

176 

193 

1634 

193 


Note that the classical Gauss elimination ends at expressions #6 and #14 and the rest must be done by 
substitution. 
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#17: nem(m, 2, 2, -1) = m 


10 0 0 


0 10 0 


0 0 10 


0 0 0 1 


BOOS 
19 B 
192S 
19 B 
176 
19 B 
1654 
19 B 


ROW_REDUCE(b) 


10 0 0 


0 10 0 


BOOS 
19 B 
192S 
19 B 
176 


0 0 10 - 


0 0 0 1 


19 B 
1654 
19 B 


This procedure offers another occasion to demonstrate that in many cases useful DERIVE- 
routines can be transferred to TI-NspireCAS. 


I 


1116 
2-113 
3 2-21 


1116 
2-113 
3 2-21 


nx. =pivot{x , 1, 1, l) 

111 6 
0 -3 -1 -9 

0 -1 -5 -17_ 


nx =pivot{nx , 2, 2, l) 

111 6 

1 

0 1- 3 

3 

-14 

0 0 -14 

3 


nx =pivot(nx, 3 , 3 , - 1 ) 

i 1 0 3 " 
0 10 2 
_0 0 1 3_ 


nx =pnvt(x i 1,1,1) 

'ill 6 
0 -3 -1 -9 

0 -1 -S -17_ 


nx =pi\--ot(nx J 2, 2, l) 

111 6 

1 

0 1— 3 

3 

-14 

0 0 -14 

3 


nx =pivot(nx , 3, 3, - 1 ) 

i 1 0 3 " 
0 10 2 
n n 1 ^ 

El 


Define pivot(v, ij, di?) = 

Func 

Local k 

For k f l,dim(v)[l] 

If dir- k<dir ■ i Then 
v(4=v[>] 

Else If k=i Then 

vM 


'pivot" stored successfully 


v 


M:=- 


v| 4 ] 

Else If dir - k>dir ■ i Then 

v[4=v[*]-M. vH 


v[v] 


End If 
End For 
v 

End Func 


David Stenenga, Honolulu 

As noted DERIVE has some nice functions doing Gauss-Jordan reduction of a matrix one step at a 
time using SWAP_ELEMENTS (v, i, j ) , SCALE_ELEMENT (v, i, s) and SUBTRACT_ELEMENT 
( v, i , j , s ) and the automatic method using ROW_REDUCE ( A, B ) . 

Since I’m teaching elementary linear algebra this semester I got thinking about this process and the 
fact that it would be interesting to see all the steps and the elementary matrices that go with each of the 
steps. As a demonstration program for my students I wrote some slight modifications of these func- 
tions and came up with the following: 
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The new functions are ROW_SWAP ( v, k) , ROW_SCALE ( v) , ROW_SUBTR ( v, k) and 
NEXT_ROW ( v) Here v is not the matrix but instead a state vector: 

[e, a, i, j , d, flag] 

where e is a vector of elementary matrices, a is the current reduction of the matrix, (i , j ) is the pivot 
position, d is used for computing the determinant and flag is set to 1 if the current column is all ze- 
ros. 

For example, you start with say [[],[0,2;2,-l ], 1,1,1 ,0] and apply the various functions using F3 or F4 
keys to put the last output into the next function. This would yield: 

(File DS-GAUSS.MTH has been preloaded as utility file.) 

#1: LOADED : \D0KUS\DNLs\DNL98\MTH32\DS-GAUSS , MTH) 




" 0 2 

i 

#2: 

[], 

_ 2 -1 _ 

H- 

O 
l 



The automatic approach is done using SHOW_GAUSS (A, B) . Here is an example: 

(\ 0 2 1 ^ 

#9: SHOW_GAUSS , I DENTITY_MATRIX(2) 

LL 2 -1 J J 

Note that the first column is just a counter, the second contains the elementary matrix which reduces 
the matrix to the one in the third column. The fourth column is the current value of the B-matrix and 
the last one is the determinant of the elementary matrix. 
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The determinant of A is just the reciprocal of the product of the fifth column entries. This is an n 3 al- 
gorithm for computing the determinant which improves the n! calculation obtained from the usual 
definition. 



/ 

'" c 

1 Z ' 



SHOW_GAUSS 



, IDENTITY_MATRIX( 2 ) 



z 

: -l . 


J 

r 

" 1 

0 " 


" 0 Z " 


" 10 " 

0 








_ 0 

1 _ 


_ Z -1 _ 


_ O 1 _ 


" 0 

1 " 


" Z -1 " 


" 01 " 

1 








_ 1 

0 _ 


_ 0 Z _ 


_ 1 0 _ 



1 

I- 1 

i 

o 


i 

l- 1 

i 


i 

o 

i 

i — l 

Z 

r-j | 


Z 


r-j | 


i 

□ 

1 

i — l 


i 

□ 

r-j 

i 


i 

l- 1 

1 

o 


# 10 : 



1 

Z 


1 


4 


1 

Z 


0 1 


" 10 " 
_ 0 1 _ 


1 1 

4 2 

1 

0 

Z 


1 


Aloha, David. 

(Another - more extended - example is included in the file DS_GAUSS.MTH, Josef) 


Albert Rich, Soft warehouse, Honolulu 

Hello Josef, 

the following might make a good DUG Newsletter item since the DERIVE on-line help does not yet 
reflect this extension of the SUB operator: 

Version 4 of DERIVE for WINDOWS (DfW) and DERIVE for DOS (DfD) introduced the SUB op- 
erator for extracting the elements of vectors and matrices. As explained in the DfW on-line help, you 
can use the SUB operator twice to extract the elements of a matrix. For example: 

[2, 3, 5; 7, 1,4; 4, 6, 8] SUB 3 SUB 2 simplifies to the 2 nd element of the 3 rd row (i.e. 6) 


" Z 3 5 " 

7 14 


= 6 


L 4 6 S J 


D 
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Beginning with version 4.06 of DfW and DfD, the SUB operator was given additional capabilities. If 
the right operand of the SUB operator was a vector of indices a vector of those elements is returned. 
For example, 

[2, 3, 5; 7, 1,4; 4, 6, 8] SUB [3,2] simplifies to a matrix consisting of the 3 rd and 2 nd row 
of the above matrix, that is 


" 235 " 
7 14 


4 6 S J 

[3, 2] 


' 4 6 S ' 

. 714 . 


Aloha, Albert. 

In the on-line help of DERIVE 6 this feature of the SUB operator is explained among the 
“Vector Manipulation Functions”, Josef. 


Parametric 3d Plots with DERIVE and 3dV 

Stefan Welke, Bonn, Germany, Spwelke@aol.com 

DERIVE does not offer 3d plots. A triplet of DERIVE functions is presented to 
produce data of parametric curve and surface plots, which can be animated 
with the freeware 3dV.exe (which is among the DNL32 files). 


Introduction 

In the DERIVE User Forum, Oscar Garcia pointed out that it is possible to animate DERIVE 3d plots 
with his freeware program 3dV.exe. I immediately downloaded 3dV and was impressed by the power 
of this tiny 57KB program and by the simple and transparent input data format of this 3d viewer. I am 
citing the 3dV.doc file: 


DATA FORMAT 

Data is given in ASCII text files, in free format, in the following 
order : 

number of points 

x, y, and z coordinates for each point (scaling is automatic) 

number of connecting lines or moves (number of items below) 
point number to draw or to move to, and color (color=0 for 
move ) 


Yes, a bit obscure, but see the sample files for examples. 

Capacity is limited by 10* (number of points) + 4* (moves and draws) 
<= 58280. 
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A triangle with vertices (1,0,0), (0,1,0) and (0,0,1) for example yields the following 3dV input file: 


3 

10 0 
0 10 
0 0 1 

4 

1 0 

2 4 

3 4 
1 4 


0 is black in 16 colour mode. We are thus jumping unseen from nowhere to the first point. The con- 
necting lines are coloured red, this is 4 in 16 colour mode. (The complete colour table is given an ex- 
pression #71 in the file T03DV.MTH.) I cannot find this obscure, this is a particular version of adja- 
cency lists, as described in DNL#26 by Hartmut Kummel. O. Garcia’s input format has the advantage 
of minimizing redundancies. By the way, 3dV is a good tool for students to design wire frames of 
solids by hand, i.e. with an ASCII editor. 


Parametric 3d plots 

It is not possible, until now, to plot parametric curves and surfaces in R 3 with DERIVE, although DE- 
RIVE apparently plots a mesh of parametric lines in 3d plots. But the user has no access to this fea- 
ture. One way around is to use projections of 3d plots to the plane. 

3dV opens another way to display 3d objects: We first produce 3d data with DERIVE , we then 
export these data to an ASCII file and finally we display our object with 3dV. This is an improvement 
to exporting DERIVE 3d plots, which are not parametric, to ACROSPIN or 3dV. 

To do the first step, we must construct functions which produce appropriate data. As DERIVE 
cannot produce data in the desired form we use lists instead. Then the triangle in the example above 
becomes: 

#1: [ [3] , [1,0,0] , [0,1,0] , [0,0,1], [4] , [1,0] , [2, 4] , [3, 4] , [1,4] ] 

Parametric plot data in this form can easily be computed with DERIVE. The following DERIVE 
code prepares the definitions of the three main functions PARAMETRIC_CURVE_DATA ( . . .), 
PARAMETRIC_SURFACE_DATA (...) and MERGE_DATA (...). 

The functions for parametric curve and surface data 

PARAMETRIC_CURVE_DATA (expr_, r_, col_) :=APPEND (APPEND ( [ [r_;3] ] , 

VECTOR (ITERATE (expr_, VARIABLES (expr_) , k_ ( r_i 2-r_| 1 ) / (r_;3-l) , 1) , 
k_, 0 , r_i 3-1 ) ) , APPEND ( [ [r_|3] ] , VECTOR ( [ j_, IF (MOD ( j_, r_| 3 ) II , col_, 0 ) ] , 
j ,l,r i 3 ) ) ) 






D 
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PARAMETRIC_CURVE_DATA ( [x (t) , y (t) , z (t) ] , [a, b, n] , col_) produces data of n 
points of the parametric curve [x(t), y(t), z(t)] for t in the interval [a, b\. This interval is subdivided into 
n - 1 parts, i.e. n equidistant points p\ = [x(a), y(a), z(a)], p„ =.[ x(b), y(b), z(b)] are calculated. col_ 
is an integer value 0 < col_ <15.0 denotes black that is no colour (back ground colour). 


PARAMETRIC CURVE DATA ( [COS ( t ) , SIN ( t ) , 0 . 3 • t ] , [ 0 , 4 tt, 1 00 ] , 15 ) 


results in a white helix (on a black background). 



As you can see, the procedure is still working - under certain conditions. My PC runs under 
WIN XP. Applying command to switch to the DOD-mode works but 3dv.exe does not work. 
Then I tried a utility DOS Box - and it works with one deficiency: The program documentation 
says that pressing any button leads back to the initial 3dV-screen. I have to end the program 
by applying the task manager. 


An example (Preload the utility file T03DV.MTH): 

#1 : "A spiral winding around a torus" 

#2 : TORUS (v) ! = [ (2+v|3COS (vjl) ) COS (v|2) , (2+v|3COS (vj.1) ) SIN (vj.2) , vj. 3SIN (v j. 1 ) ] 
#3: LINE (t,m, r) := [t,m-t, r] 

#4: PARAMETRICCURVE DATA (TORUS (LINE (t, 1/12, 6/5) ) , [0,24 tt, 400] ,15) 
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P ARAME TRI C_S URFACE_DAT A (expr_, xr_, yr_, col_) : = 

APPEND (POINTS_ (expr_, xr_, yr_) , LINES_ (xr_, yr_, col_) ) 

PARAMETRIC_SURFACE_DATA ( [x(s,t) ,y(s,t) ,z (s,t) ] , [a,b,n] , [c,d,m] , [cs,ct] ) pro- 
duces data of the parametric surface \x(s,t), y(s,t), z(s,tj] with 5 e [a, b] and t e [c,d\, with n points in 
interval [a,b] and m points in interval [c,d\. cs and ct specify the colours of the 5- and t-parameter 
curves respectively, so the two parameters can be distinguished graphically. If you don’t want differ- 
ent colours for the parameter lines, then you can use [c] as fourth argument. 

The order of the intervals is prescribed by the built-in function VARIABLES. For example 
VARIABLES ( [s, s*t / '2, t] returns [s, t] , so the first interval belongs to the parameter s, the 
second one to t. 

The next example uses two colours: 

#5: PARAMETRIC_SURFACE_DATA (TORUS ( [s,t,l] ) , [0,271,12] , [0,271, 30] , [8,15] ) 

The figures displaying the spiral and the to- 
rus are original figures from DNL#32 (1998). 

All other figures are from 2014. 


At the end of Stefan’s contribution we will 
generate the objects with DERIVE 6 and 
Tl-Nspire as well (page 15). 

Combining data of two objects 

MERGE_DATA ( v_, w _) : =F2_ ( v_, w_) +F3 (w_, v_) 

MERGE_DATA (datal , data2 ) combines data of two distinct pre-3dV objects to one single 
pre-3dV data object. Now we merge the torus and the spiral: 

#6: MERGE DATA(#5,#4) 




This picture is made 2014 applying again the DOS Box (page 21) and the DERIVE file from 
1998. I followed Stefan’s instructions and produced a screen shot. Finally I reversed the col- 
ours using a graphics program, Josef. 
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Procedure of conversion to 3d-format 

PARAME TRI C_CURVE_DATA , PARAMETRIC_SURFACE_DATA and MERGE_DATA produce 
vectors of a structure which is appropriate for conversion to 3d-format of Garcia’s 3dV (= 3d Viewer). 
In order to convert the data for 3dV proceed as follows: 

You first need the free awk-interpreter mawk.exe which is described by O. Garcia in his DNL#25 
contribution. Then you need the script file convert . awk which is a simple ASCII-file: 

convert . awk 

# mawk -f convert. awk datei.in > datei.out 

{RS="\] z=NF ; 

{ for (i=l ; i<=z ; i++) 

{ 

getline < FILENAME 
gsub ( "\*POW\ (10, ", "E") ; 

gsub ( ,f \) ", "") 
gsub $0) 

gsub (","," ",$0) 
print ( $0 ) 

} 

} 

} 

convert . bat 

REM Conversion batch file infile. pas -> outfile.3d 
mawk -F "\] n -f convert. awk %1 > %2 
del %1 

After DERIVE has calculated the data, save the resulting data- vector as a PASCAL-file (File > 
Write > Pascal file ...), say data . pas. Then call the batch file convert . bat with the input file 
data . pas and say figure . 3d as arguments: convert data, pas figure 3d. Now you are 
ready to load figure . 3d with 3dV . exe and animate the figure with the mouse. 

mawk . exe, convert . awk and convert . bat must be in the same directory (folder). Notice 
that the original *. pas-file will be deleted after conversion. 

Limitations 

Limitations are due to three reasons: 

1 Memory of DERIVE. The DOS version does better that DfW. Too many points produce too long 
lines in DfW. In the DOS version the limit on my machine seems to be a 15x 15 grid for the sur- 
faces and 300 points for curves. 

2 Memory usage of the mawk interpreter when converting the *.pas files to *.3d files. 

3 The maximum number of 10*(number of points)+4* (moves and draws) < 58280 for 3dV.exe. 

In my practise the first item is the real restriction while the last one is only a theoretical restriction. 
I assume that DERIVE XM overcomes these difficulties. Because I have no access to this version 
I calculated the data for the examples in this paper with Mathematica. 


MAWK strips brackets and commas 
off in the input . pas file. Finally you 
need a batch file convert . bat which 
manages the conversion job: 
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(/ tried with DERIVE XM and could reproduce Stefan’s examples without any problems. 
From earlier experiences with similar problems (ACD, ACROSPIN) I recommend to reduce 
the DERIVE-output to 3 decimal digits. DERIVE versions 4 for DOS and WINDOWS as well 
didn’t accept the MERGE command. Trying to overcome this obstacle I found out that in two 
auxiliary functions several times a SUB[1,1] was used which worked in earlier DERIVE ver- 
sions. Now you have to write SUB 1 SUB 1 instead. I sent a message to Stefan before in- 
specting my mailbox and then found his message that he has found out that ..., so we had 
the same idea nearly at the same time. It might be of interest for you to change your DE- 
RIVE code if you want to use well-tried DERIVE files with an upgraded DERIVE version. 
See also Albert Rich’s message in the User Forum. Josef) 


Further examples 

I close with two examples. Dini’s surface is a twisted pseudosphere. This is a surface with constant 
negative curvature. Its explicit parameterization is: 

DINI (s, t) ; = [COS (s) SIN (t) , SIN (s) SIN (t) , COS (t) +LOG (TAN (t/2) ) +0 . 2 • s] 
PARAMETRIC_SURFACE_DATA (DINI (s, t) , [0, 471, 30] , [0.01,2,20], [4,1]) 

The sea shell is defined as follows: 

SHELL (r, s) : = a-r- [COS (r) (b+COS (s) ) , SIN (r) (b+COS (s) ) , c+SIN (s) ] 

[a «= 1, b >= 1, c >= 1.3] 

PARAMETRIC SURFACE DATA ( SHELL ( r, s ),[ 0 , 2n, 30 ] , [ 0 , 371, 30 ] , [ 4 , 5 ] ) 



Originals from 1998 
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2. 3DV25.zip and mawkl22.zip can now (2014) be found at 
http://cd.textfiles.com/simtel/simtel20/MSDOS/GRAPHICS/ 
http://www.filewatcher.eom/m/mawkl22x.zip.102254-0.html 

3. A. Gray, Moder Differential Geometry of Curves and Surfaces, CRC PRESS, 1993 


and among the DNL#32 files 
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DERIVE 6 (above) and TI-NspireCAS (below) 
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AN EXERCISE ON THE EVALUATION OF THE MODIFIED 
BESSEL FUNCTIONS IN DERIVE 

S. Abraham Ibrahim and P. Fernandez de Cordoba 
Departamento de Matematica Aplicada,Universidad Politecnica de Valencia, Spain 

Abstract 

Using the continued fraction method in DERIVE we present an algorithm to calculate 
Modified Bessel functions of fractional and integer order based on this method. Us- 
ing Derive we also compare with the traditional Miller s algorithm. DERIVE program 
files and some illustrations of the results are included. 

1. Introduction 

When a given function is approximated using the continued fraction representation, we observe that it 
generally converges faster than other approximations i.e. Taylor expansion. We insist on this in sec- 
tion 2, where we also present as an example the continued fraction representation of the exponential 
integrals. In sections 3 and 4, we present the algorithm to calculate Modified Bessel functions (MBFs) 
of integer and fractional orders, based on the continued fraction method. 

We compare our algorithm (which as we show is very efficient), with the usual numerical methods (in 
section 5) which calculate MBFs taking into account normalization relations (1,2) . 

In this paper we introduce an algorithm and corresponding DERIVE codes to evaluate regular and 
irregular MBFs without any re-calculation through normalization relations. Furthermore, the method 
maintains the stability of each recurrence relation, i.e., we use forward recurrence relations for the 
MBFs of the second kind and backward ones (3) for the MBFs of the first kind. The algorithm uses for- 
ward recurrence relations to generate irregular MBFs and takes into account the continued fraction 
method to evaluate high order regular MBFs. From these values we can generate regular MBFs apply- 
ing backward recurrence relations. 

In Section 6 we present the structure of the programs. A comparative table with some results is pre- 
sented in section 7. 

For the importance of DERIVE software as support for teaching and researching in science, we evalu- 
ate the continued fractions representation of a function and its application to the calculation of MBFs 
with a code to be used with DERIVE. 

2. Evaluation of continued fractions 

A continued fraction is given by (4,5) 




( 2 ) 
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We introduce the SMBFs of the third kind, k n (z ), and the first kinds, i n (z) as particular solutions of the 
differential equation 

z 2 co\z) + 2 z co\z) - [z 2 + n(n + 1] a>(z) = 0 

In the code, wie calculate SMBFs of all orders below N max , i.e., we generate the set 

SMB{z) = { i n (z ), k n (z); n = 0, 1, 2, ..., N max ) 

• Evaluate all the SMBFs of the third kind, { k n {z ), n = 0, 1,2, ..., N max }, taking into account the 
known values of k 0 (z) and k\(z) using the forward recurrence relation: 

k n+1 (z) = ^-^-k n (z) + k n _ l (z). (13) 

z 

• Use the continued fraction method to evaluate the ratio 


/ ,(*) 


FC(z) = - 


V' (z) <z> 2 r« + — ) 
2 l 2j 


2 + - 


2 w + : 


• Calculate the upper order SMBFs of the first kind, i Nmax using the already known values of 
k N (z) and k N (z), the ratio FC(z), and the value of the Wronskian of SMBFs (5) 

W ^ k ^-, = ( z ) ( z ) + ( z ) {z') = nt 2 z 2 

As in the previous expressions, we can write: 


and then, 


h 0) = / 2z 2 [FC (z) ^ (z) + (z)] 


iv , ( z ) = FC(z)i N (z). 


Notice that we have calculated not only i N (z) but also / iV (z). 


Evaluate all SMBFs of the first kind, {4(z), = 0„ 1 ,2, ..., N max }. taking into account the calculated 

values of z /Vwax (z) and i fc i(z) using the backwards recurrence relation 

i n _ l (z) = ^-^i n (z) + i n+l (z). (18) 

z 


We would like to point out that / 0 (z) and fi(z) can be used as checks on the accuracy of the procedure, 
because both are calculated at the end of the method. 
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5. Miller’s algorithm 

Miller’s algorithm initiates' 1 calculating directly an sufficient high order m SMBFs of the first kind 
(m » N max ) with the values of “0” and “1” respectively for i" m (z) and z"_i (z) using the backward re- 
currence relation 

i’n - 1 ( Z ) = — O) + Cl ( Z ) ( 19 ) 

Z 

i.e., we generate the set \i rf n (z), ?2 = 0,1,2, ...,mj. 


Taking into account the calculated values of /"(z) and /" (z) and using the known values 3 * (5) 
and /, (z) we calculate the re-normalization constant 


i',(z ) 

Finally we generate the set 


i n (z) = p(z) (n = 0, 1 , 2, ... A«ax) 


of z 0 (z) 


( 20 ) 

( 21 ) 


6. Program specification 

In this section, we introduce three DERIVE programs. A first code to evaluate the SMBFs using 
Miller’s algorithm, a second one to calculate the continued fractions representation of the integral 
function and the last one to evaluate the SMBFs using the continued fractions. We pass to detail the 
structure of these programs. 


1. DERIVE program to evaluate SMBFs using Miller’s algorithm 

(see appendix A) 

The code is organized as follows: 

- with expressions #1 and #2 we calculate the known (5) values of /o(z) and z'i(z), 

- with expression #3 we generate the set i” r , (n = 0, 1,2,..., N max ), 

- with expression #4 we evaluate the re-normalization constant, 

- and finally with expression #5 we generate the set /" , (n = 0, 1,2,..., N max ) 

2. DERIVE program to evaluate the continued fraction representation of the integral 
function (see appendix B) 

- with expressions #1 and #2 define the a’s and b’s in expression (4), 

we evaluate the values of a n and b n respectively with the functions A(k,n) and B(k,z ), 

- define the d s and yS’s in expression (5) using the recurrence relation (6) and (7), 
we evaluate the values of cti and fij with the functions a(k,n,z) and fi(k,n,z ), 

- we store the values of fk = ccklfik , (k = 0, 1,2, ..., N) in the H(£,?/,z) vector 

3. DERIVE program to evaluate SMBFs using Miller’s algorithm 

(see appendix C) 

- with expressions #1 to #1 1 we calculate the initial values of K 0 (z) and K\(z) using the poly- 
nomial approximation (5) . For this we use the functions AQ(z) and K\{z). 

- with expression #13 we evaluate the set { K n {z ), n = 0, 1, 2, ..., Nmax} using equation (8) 
with the functions FBK(n,z). 
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- with the expressions #14 to #18 we evaluate the continued fraction storing it in the 
//(#,?z,z)-vector. Finally we calculate it using the FC(k,n,z)- function. 

- Evaluate the set {4(z), « = 0, 1,2, ..., Nmax} using function FBI(k,n,z) in expression #19. 


Unfortunately I did not receive the respective files for 2. and 3. in 1998. I wrote a mail to the 
authors in 2014. They answered but they could not retrieve these files, Josef. 


7. Results 

In this section, we have presented some numerical examples of the accuracy of the method. 

1 . Results of the program to evaluate the integral function E n (z ) using a continued fraction repre- 

sentation. 

We present the results for some values of “n” and “z” in the following table using the fraction 
continued approximation: 


Table 1. Calculated values of the integral function E n (z) using the continued fraction method. 


k 

n 

z 

E n (z) 

Abramowitz and Stegun 

20 

20 

0.10 

0.047359994 

0.0473600 

20 

10 

0.69 

0.052505381 

0.0525055 

25 

4 

1.15 

0.071164022 

0.0711632 


2. Results of the evaluation the MBFs are represented in the next table:. 

Table 2. We compare our results with those implemented in DERIVE’ s files. 


n 

z 

MILLER 
(1=5, m=5) 

CONT. FRACT. 
Method 

k 

BESSEL_I(ft,z) 

BESSEL I 
SERIES (ft ,z,/ft) 
(m= 10) 

Abramowitz and 
Stegun 

1 

5 

24.335642097 

24.3356421099 

10 

24.3356421424 

24.335642065 

24.33564214 

5 

5 

2.157973762 

2.157974553 

10 

2.157974547 

2.157974547 

2.157974547 

10 

2 

3.016963885*10"' 

3.016963839* 10' 7 

11 

dub. Accuracy 

3.016963878* 10' 7 

3.016963879*10‘ 7 

15 

2 

8.139432548*10' lj 

8.139432410*10 13 

5 

dub. Accuracy 

8. 1 39432530* 1 0"' 3 

8.139432531 *10“ 13 


8. Conclusions 

We have introduced the DERIVE version of an algorithm to evaluate the continued fraction representa- 
tion of a rational approximation of a function and an efficient algorithm to generate MBFs. This code 
illustrates the se of the continued fraction method and gives an accurate way to evaluate the MBFs 
although it is slower than the implemented BESSEL. MTH file. 
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Logos, Stars and other Figures with 
DERIVE and the 77-92 

Josef Bohm, Wurmla 

I am a teacher at a Business College and our pupils come to our school when they 
have finished lower secondary education at different types of schools. So one of my 
tasks is to find a common base to teach mathematics repeating basic mathematical 
knowledge. As we have to deal with functions and their graphic representation it is 
obvious to freshen up the system of coordinates. Some years ago I did this success- 
fully supported by DERIVE and now class room bounded I am using the TI-92. Us- 
ing the DATA-tool simultaneously with the Data plot on the Graph screen it is very 
easy and entertaining to explore the plane divided by the axes. 

I could see the students find transformations like translations and reflections by 
themselves. The points can be connected to draw figures. 

Draw a square which has one vertex in each quadrant by defining its 
points. (It is always nice to see them explore how to close the polygon) 

Draw a house. Using more pairs of data columns you can show more sepa- 
rated objects. 

Draw an equilateral triangle with a given base line. (They have to remem- 
ber the formula for the altitude and work exactly using V3) 

Complete the triangle to a regular hexagon. 

Shift the hexagon to another midpoint. 

(Some students had the idea to apply the translation(s) immediately on the 
data columns: c3=cl + 2 and c4=c2 - 4 produces a translation by the vec- 
tor (2,-4)). 

If you are interested in the work sheets then please write, email or call. The paper is 
in German, but it is very easy to follow even in other languages because there are a 
lot of pictures. If there are many of you interested I would try to produce a transla- 
tion. 

See two short parts of the work sheets: 

4. Wechsle mit wieder ins Datenblatt und 

iibertrage nun Punkt fur Punkt die Punkte A, B, C 
und D aus Aufgabe 2 ins Grafikfenster. 

Switch to the DATA-MATRIX-Editor and plot points 
A, B, C, and D from task 2 into the GRAPH-Window. 
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1 1 . Zeichne den Umriss eines Hauses im 1 . Quadrant. 

(Plot the graph of a house into the 1 st quadrant.) 

(<- nur ein mogliches Beispiel) 

Uberlege zuerst die optimale Reihenfolge der Punkte im 
Polygonzug. 

Consider the best order of the points for plotting the poly- 
gon. 

12. Spiegle dieses Haus zuerst an der x-, dann an der y-Achse und abschliefJend am Ko- 
ordinatenursprung. Gib fur jede Spiegelung die Konstruktionsvorschrift an. 

(Plot the mirror image of your house wrt the x- and to the y-axis and then wrt the origin. 
Describe the respective construction rule.) 

Most of this can be done with DERIVE and the TI-92 as well. 

At last I gave a home work: 

The students should design their own personal logo, their initials, a fantasy 
star or any other object they liked to represent in the plane using coordinates 
to describe it point for point. 

I promised to provide a program for their devices which would animate their crea- 
tions. 

20. Als Abschluss sollst Du Dir Dein eigenes LOGO - Initialen oder sonst etwas - im Koor- 
dinatensystem schaffen. Erzeuge ein neues Datenblatt LOGO. Jedes Polygon (Vieleck) 
wird durch ein Spaltenpaar beschrieben. 

Als Muster siehst Du hier ein “E.”. Die Ach- 
sen, ihre Bezeichnung und das Gitter wer- 
den am Ende ausgeblendet. Du kannst da- 
mit rechnen, den ganzen Schirm zur Verfii- 
gung zu haben. 

(-11 <x< 11, -5<y<5) 

Wenn das LOGO fertig ist, werden wir es mit Hilfe eines Programms animieren, d.h. einen 
kleinen Film draus machen. 

Mr IN LOGO 



MiilH 




Fiffp E ii H ■= T 


FUN-1 
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We - my colleague Tania Koller and I - gave this homework in four classes and we 
didn’t expect too many nice results. But sometimes even teachers have their won- 
derful experiences. After the week end when we asked the pupils to present their 
designs using the view screen our eyes became bigger and bigger: some of the boys 
and girls might have spent hours to realize their ideas. Next day I took my notebook 
with me and using GraphLink collected the Four Classes Picture Gallery. You can 
see the Ha collection below. 





r i_n 

cl _ r 

LnJ 

LnJ 


And nobody asked: why are we doing that in such a “boring” way? We could use 
Paintbrush, Corel Draw or any other painting program. They had just fun and 
proudly presented their screens. Surprisingly some of the weak mathematicians 
had excellent ideas and could realize them. 












D 


-N-L#32 


Josef Bohm: Logos, Stars and other Figures 


p 35 


But then they reminded me on my promise to make a “movie” of their pictures. 


In pairs of columns I produced letter after letter according to a design on a grid pa- 
per and wrote the data into a data table called "1 ogo". Then using F2 Plot Setup 
I called FI Define and there defined one polygon after the other. The screen shot 
shows the creation of the data plot of columns 7 and 8 representing the contour of 
"9". 12 columns make the whole picture. 



\*m 


imii 


DATA 









c 1 

c 2 

g3 

G 4 

gS 

g6 


1 

“8 

*4 

“3 

-4 

■i icgn; 

0 


2 

*8 


-3 

4 




3 

-6 


tns.-. 

■4 

PMol 



4 


4 

EMU 


BM 

mm 

■ 

Fb 

BIVU 

4 

*3 

-4 

Bgj^PTPl 

□ 

■ 

£ 

mm 






■ 

7 

ieaiMfS 






■ 


cl .Iitle= 

i rnn ' Rnn niun fuhc a 


I 



Plot Type- ....... 

xyl ine* 


Hark, ............ 

Dot 4 


x 

c? 

*s 

if 

y — - - 

i-il-sL- Z<Cz-yxt-K Mjdtr; 

p 


Use Preq and Categc 

ries? HQ-» 


r ^ ;i’. try r .» [i* s . ^ ^ ^ ■ 



Cr=c i4 (rtw x*hi- 

2- 




if i n i iii {"ftj 


,m ; E » [EN7Eft3-DK ft HI) rF^~l=rftHi~n 



I- • •• - ■■ — ~ — ■- 1 Mil 

1 

1 

a 



p l;f 






Liii m *wth riiNca : 


I wanted to produce a movie merging the 
TVs logo and "DERIVE. So I had another 
data plot with the data fixed in variable 
"logod". 

Now I needed a program BILDER( ) to 
produce the pictures which should be 
"cycle-pictured" by the next program to 
perform the animation: 

I came across some problems: I was not able to address the variable "1 ogo" - the 
name of the data table - as a program parameter. I wanted to call BILDER(1 ogo) to 
make the program more flexible - also for BILDER( logod) - but I failed. The same 
reason made it impossible to use dim (logo) etc, but fortunately and inconse- 
quently 1 ogo [ i ] returns the i-th column of the data. So I used the number of col- 
umns as parameter. The program produces three families of pictures "dr"s, "ki "s 
and "st"s. You see the pictures appearing on the screen, each one is saved in one 
extra file. For my students I made a version including the choice how to represent 
the movie: 


( 

i 

) 

Rotate 

Rotation around the vertical axis 

( 

2 

) 

Swi tch 

Rotation around the horizontal axis 

( 

3 

) 

Flash 

Zooming in and out 


The second program f i 1 m ( ) does the animation. 
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The smaller versions are on the diskette (BILDER1, FILM!). The super version for 
merging the TI-92 and the DERIVE-Logo is also on the diskette (BILDER2, FILM2). 

You can work very similar using DERIVE. I have three different ways to show the 
figures (see below). Try and be surprised: logo*.acd-files and logo.mth. You have to 
use ACROSPIN in combination with ACD (Diskette of the year 97, DNL#28). In- 
spired by Stefan Welke’s 3dV-contribution I produced a DERIVE-function to trans- 
fer polyhedrons into 3d-format. You can find the logo*.3d-files on the diskette. 

Bi 1 der{ 6n ) 

Prgm 

Local i , j,et,xl ,yl 
Cl rGr apli ; FnOf f 
ZoomDec 

setGr aph ( "Axes” , "Of f M ) 
setGr a ph ( "Grid”, “Off" ) 
setGraphC "Labels" , "Off”) 

For 0t , 0 
Cl rDraw 
For i , 1 , en , 2. 

1 ogo[ i ]->xl : 1 ogo [ 1+1 ]-»yl 
For j , 1 , dimtxl ) -1 

Line xl [ j]*cos(et) ,yl [ j] ,xl [ j+l]*co$ 

EndFor 

EndFor 

StoPic #( "dr , ’&stri ng( 10*et/(*)+l) ) 

EndFor 

For 

Cl rDraw 
For 1 , 1 ,ert, 2 

1 ogo[i ]->xl : 1 ogo[ i+1 ]+yl 
For j , 1 , dimtxl )-l 

Line xl [J] ,yl [ j]+costet) ,xl [ j+l] t yl [ j+l]*cos(et) 

EndFor 

EndFor 

StoPic #( ”ki ” &$tri ng ( 10*et/ (x )+l ) ) ■ 

EndFor 
For flt,l ,11 
Cl rDraw 
For 1 , 1 , en ,Z 

1 ogofi 1 : 1 ogo[ i+1 ]-»yl 
For j . 1 , dimtxl )-l 

Line xl [ j ]*< 11 -9t )/ 10 ,yl [ j ]*( 11 -et ) / 10 ,xl [ j+l]*( 1 1 -et ) /10 , y 1 [j+l]*Ul- 

eti/io 

EndFor 

EndFor 

StoPic # ( " st "&st M ng ( 8t ) ) 

EndFor 

EndPrgm 




filmO 

Prgm 

Lbl start 

CyclePic "st" , 11 , . 12.1 , -1 
CyclePic " drM 1 , . 12 , 1 , -1 
CyclePic "111 ",11 r . 12, 1,-1 
Goto start 
EndPrgm 


(et).yl [J+l] 
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The DERIVE - graphs are based on the logo given below. The file is among the ac- 
companying files of the DNL. 


logo 


-2 -3 " 

-2 -2 
-2 3 

-6 3 

-6 5 

0 5 

0 -4 

-1 -5 

-5 -5 

-6 -A 
-6 -2 


" 0 5 " 

5 5 

6 4 
6 1 

5 0 

6 -1 

6 -4 

5 -5 

0 -5 


' 23 " 
4 3 

4 1 

2 1 
. 23 . 


" 2 -3 " 

4 -3 

4 -1 

2 -1 
_ 2 -3 _ 


-4-2 L 0 5 

-4 -3 


-2 -3 



More interesting for the revised version of DNL#32 is how to perform the movie on 
TI-NspireCAS. The coordinates of the logo points are given in lists in a Lists & 
Spreadsheet page. Take care that columns are of equal length, i.e. fill up empty 
cells with the coordinates of the last point (see II 1 and 12 2) . 

Collecting all lists in another list (in the Calculator) results in a matrix (see next 
page). I called this matrix logo. Three small programs produce list II, 12, ..., In with n 
= number of columns of the matrix (= 2 x number of polygons). 



AH _1 

B 11_2 

C 12_1 

□ 12 2 

E 13_1 

F 13 2 

G 14_1 

H 14 2 

| 

- 










1 

-2 

-3 

0 

5 

2 

3 

2 

-3 


2 

■2 

-2 

5 

5 

4 

3 

4 

-3 


3 

-2 

3 

6 

4 

4 

1 

4 

-1 


4 

■6 

3 

6 

1 

2 

1 

2 

-1 


5 

■6 

5 

5 

0 

2 

3 

2 

-3 


6 

0 

5 

6 

-1 

2 

3 

2 

-3 


7 

0 

-4 

6 

■4 

2 

3 

2 

-3 


8 

-1 

-5 

5 

-5 

2 

3 

2 

-3 


9 

-5 

-5 

0 

-5 

2 

3 

2 

-3 


10 

-6 

-4 

0 

5 

2 

3 

2 

-3 


11 

■6 

-2 

0 

5 

2 

3 

2 

-3 


12 

-4 

-2 

0 

5 

2 

3 

2 

-3 


13 

-4 

■3 

0 

5 

2 

3 

2 

■3 


14 

-2 

-3 

0 

5 

2 

3 

2 

-3 

1 

E 

IE 


Cn.D14\ 
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logo ={ 

n in 

2,12 

1,12 

2,13 

1,13 2,14 1,14 2 } 




0 

rot_x 

6 16 


-2 -2 

-2 

-6 

-6 

0 

0 -1 -S -6 -6 

-4 

-4 

-2 


Define rot_x(p?c)= 

1 

B 

n 


-3 -2 

3 

3 

S 

5 

-4 -S -S -4 -2 

-2 

-3 

-3 


Prgm 




0 S 

6 

6 

5 

6 

6 5 0 0 0 

0 

0 

0 


Local d,i 




5 S 

4 

1 

0 

-1 

-4 -5-5 5 5 

5 

5 

5 


d =dimjpic)[l] 




2 4 

4 

2 

2 

2 

2 2 2 2 2 

2 

2 

2 


For i, l,d,2 




3 3 

1 

1 

3 

3 

3 3 3 3 3 

3 

3 

3 


mat ► listjpic [? ])-»#( 1 1 1 " &strin g (?) ) 




2 4 

4 

2 

2 

2 

2 2 2 2 2 

2 

2 

2 


matHist(p?c[?+l]- cos (*■ 2- ji)) -^#("l"&string(?+l)) 




-3 -3 

-1 

-1 

-3 

-3 

-3 -3 -3 -3 -3 

-3 

-3 

-3_ 


End For 



rot_x(logo ) 







Fertig 

End Prgm 



DelVar t 







Fertig 




rotjy{logo) 







Fertig 




DelVar t 







Fertig 




_ . A 










0 



B 


rot_y 0/6 

"flash" erfolg. gespeichert | 

Define rot_y(p?c)= 

Prgm 
Local d,i 
d. =dim^?r)[l] 

For i, l,d,2 

mat ► lister [?]■ cqs(l 2- ji)) -*#( "l"&sting(?)) 
matMi5t|pjc[i+ 1]) -* #( "1 " &string(?+ 1)) 

End For 
End Prgm 

Define flash^??c)= 

Prgm 
Local d,i 
d =dim^?r)[l] 

For i, l,d,2 

mat^list(p?c[?]- sin(f- 2- ji)) -*#("l"&string(i)) 
matMist(pjc[i+l]- sin(/- 2- ji)) -»#("l"&string(?+l)) 

End For 

End Prgm U 


After running one of the three programs insert a Graphs & Geometry page, intro- 
duce a slider for parameter f with 0 < t< 1, step size 0.01, and enter the lists II 
through In in the Scatter Plot entry lines. Moving or animating the slider shows the 
“movie”. 

If you want to have another movie, first delete t in the calculator, run the next pro- 
gram, switch to the Geometry page and initiate again the slider. Have fun! 
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J J 

3.22 y = ;.sin$ V + — ,[oosi|>.V + sin$. V] + — ,[2.cos$.V + sinf.fv - v)j + 


/4 j (d 4 V Y1 

+ - cos 4-(5V - V) + sinf \^J - 3Vjj + 


t 5 '' ( d 4 V j,'' 

»r“*W ' V J + 


■w 

{ 5^ 


r^v ... .v 

.1 t - d V + V I + 


r^v ,„^v 


+ «■ + V J + sin *W ■ w i ? + 5 T] + 



f 

1 

f>V 


.3 

■ J 

d 7 \ 

,-d S V „« 

-+- 

7 f 

COSO'. 

6 x 

V rf<|> S 

- 20 — — 
</( i>* 

+ <> vj + 


d $ 7 

- 15 - c + 15 V 


t s \ 


f d 7 v 

d s V 


•l 1 



+ 


COSO. 

7 — 

- 35 ~nj 

+ 21 V - 

V + 




8 ! \ 

V <i<|/ 

d$ i J 


J] 





- 

^/ s v 

t/^V 


A] 



+ 


sm$. 

0 

2 / — 7 - + 

55 , ■ 

■ ?v 

+ 



8 !' 



d § 4 

)\ 




t 9 f d S V d 6 \ d 4 \ 

— .cosf 5 j - 56— J + 56— J - S\ i 

91 dty 8 dif 6 dif 4 ) 



d 9 V 

-,J 7 V 

, 

w* sm * 

U<t> p 

28 7 + 

d$ 7 


(d 9 \ 

d 7 V 

——.COSO 

9—jr 

- 84 _ 

10! 1 


rf (|> 7 

ll ° ■ , 

( d I0 \ 

3 je rf*V 

lor sm * 

\di IB 

rf«i» s 


126 — -r - 84— r + PV 
d§ 6 d§ 4 


Whewell Type 2, ® = W(s) = W J5 requires 

dx d 2 x_ ci^jc d n x ^ ^ dy_ d 2 y d J y d n y 

ds‘ ds 2> ds 3 ds n ds’ ds 2 ds 3 ds n 


To use DERIVE declare W = W(s), differentiate with respect to 6* using DERIVE. 


4.00 <t> = W (s) = W 4 = w 
dx 

4.01 — ~ cqs<|> = cosW 
ds 


d 2 x 

4.02 — r = -W'.sinW 
ds“ 


4.03 — T - - W".sinW - (W') 2 .cosW 
ds * 

• ^ 

4.04 — j - -cosW.fW' + 2WW'J + sinW.[(W ') 2 - W*] 


4.11 — = sin<t> = sinW 
ds 

d 2 y 

4.12 — £ = W'.cosW 
ds 2 


4.13 —7 = W'.cosW - (W') 2 .$inW 

d s* 


4.14 - -cosW.[(W ') 2 - W'l - sinW.fW" + 2WWJ 

d v 4 1 


/i f* 

4.21 x = t - — .(JP ) 2 - — .[W' + 2W'W"] 

6 24 

*2 * 3 *4 

4.22 y = — W' + — W” - “ |(Wr - W w 


where (f) = 0 
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I 


Euler Type 1 is p = E(O) = E<d, requiring 
dx d 2 x d 3 x d n x ^ d 2 £ d 3 y_ d n y 

d-V d§ 2> rfij) 5 d$ n d§ 2 ' dfy s d§ n 


To use DERIVE declare E = E(O), differentiate with respect to <D using DERIVE. 
This leads to: 


5.21 ^ 


R t + 



r| 

6 


(E - 3b).t 4 

24 


5.22 y = 


E .t 2 E.t 3 

+ + 


(3 E - E).t* 

24 


where 0 = 0. 


Jl 


Euler Type 2 is O = F(p) = Fp, requiring 

dx d 2 x d 3 x d n x , dy d 2 y d 3 y d n y 

— — - 7 ... and — — t, ... — — . 

dp d p 2 dp 3 d p” dp dp 2 dp 3 d p n 


To use DERIVE declare F = F(p), differentiate with respect to p using DERIVE to obtain: 


6.21 * = /.| 

6.22 y = - 
where 0 = 0. 


r rfFi ' 

(P) 

d F d 2 F 


~d 2 F d 3 F 

d F 

i p >; + 

ui 

.dp P d p 2 . 

+ 

1 Ui 

2. . • p. , 

. dp 2 dp 3 

P dp. 


r dT~ 

(EvJ 


>3} 


d F 


2 A dp 


d 2 F 

"’Vj 


TABLE OF OSCULANTS 


(1) Quality, (2) Order. (3) Number, (4) Equation Identity 


a) 


(3) 

NAMES(S) 

EQUATION 

0) 

Qi 

1 

1 

Gradient 

— - tan^ 
dx 

Ul 

QU 

1 

2 

Elasticity 

_y^_ 

x dy 

6.10 

Q2 

2 

I 

Radius of Curvature 

dx 

2.20 

Q2a 

2 

2 

Gradient Change per 

Abscissa 

increase 

__ secV 
dx 1 p 


Q3 

3 

i 

Spiral ation 

dp p 

-E - p y = tans = — 
ds p 

4.11 

Q3a 

3 

2 

Aberrancy or 
Deviation 

i d P = 4 = £ = KnS 
3 ds 3 3p 

4.10 

Q3b 

3 

3 

Involation 

dp 

T* =p = #> 

19.0 
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(1) 

(2) 

(3) 

NAMES(S) 

EQUATION 

(4) 

<?3c 

3 

4 

Whewell 

Gudermannian 


16.6 

Q3d 

3 

5 

Spiralation-2 

{*_ = P_ 

P P 1 

20.0 

Q3e 

3 

6 

Clothoidation 

P' ^ P 

P 1 " 

21.0 

Q4 

4 

1 

Radius of 
Tau-ature 
Radius of Q4 

__*_/ + (p'f„ _ p(p + P 2 ) 

da p pp - p 2 

22,0 

Q4a 

4 

2 

Catenation 

d 2 P _ _ PP - 

di P p 3 

9.10 

Q4b 

4 

3 


d p - p - pj(p') + "*'} 


Q4c 

4 

4 

Cycloidal 

Displacement 

l=pp‘ + (p’f 
P 

23,0 

B 

4 

5 

Tractration 

\ 

ii 

24.0 

Q5 

5 

i 

Betation 

* (p") 

27,0 

Q5a 

5 

2 


_ P 4 P - 4p 3 pp + 1P 2 P 3 
P 


Q6 

6 

1 

Radius of Q6 

ds 

Tp 


Q7 

7 

1 

Gammation 

i 

i 


Q8 

8 

1 

Radius of Q8 

ds 

dy 







There are different scalings of the 4 th order Curve of constant Tau-ature, which is the forth in the ho- 
mologous series of curves of constant geometric quality, whose first three are the Straight Line, the 
Circle and the Logarithmic Spiral. 
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Leibnizens silbeme Taschenuhr 
Leibniz’ Pocket Watch of Silver 

Josef Lechner, Viehdorf, Austria 


In H. Heuser’s worth reading book “Gewohnliche Differentialgleichungen” [1] you can find the fol- 
lowing problem as an introduction (Problem of the Tractrix - Drag Curve, Curve of Pursuit): 


Leibnizens silberne Taschenuhr oder die Traktrix 

Gottfried Wilhelm Leibniz (1646-1716) behandelte 
1693 in den Leipziger Acta eruditorum das folgende 
Problem (s. Abb.): in der xy-Ebene ziehe man einen 
Punkt P an einer straff gespannten Schnur PZ der 
Lange a. Der Zugpunkt Z soil auf der y-Achse fort- 
riicken, und zu Beginn des Vorgangs beflnde sich P 
in ( a,0 ) (= „Anfangsbedingung“). 

Welche Kurve beschreibt P? 

Des „Verstandnisses wegen” imaginierte Leibniz eine horologio portatili suae thecae argenteae , die 
an ihrem Kettchen per tabulam gezogen wird. Erfinder des Problems sei Claudius Perralus, Medicus 
Parisinus, turn et Mechanicis atque Architectonis studiis egregius ... Gemeint ist der Pariser Architekt 
Claude Perrault (1613-1688), der die beruhmte Saulenfassade an der Ostfront des Louvre entworfen 
hat; die Heilkiinste des Mannes diirfen wir auf sich beruhen lassen. 

Da die Zugschnur PZ zur gesuchten „Zugkurve“ oder „Traktrix“ y = y(x) offenbar immer tangential 
ist, konnen wir aus der Figur sofort die Gleichung 

, si a 2 - x 2 

y = 

X 



ablesen. Da in ihr die Ableitung y der gesuchten Funktion auftritt, ist sie eine Differentialgleichung, 
wenn auch nur von de extrem einfachen Form y = /(x). Bereits die elementare Integralrechnung 
lehrt, dass alle Funktionen 



n 2 

si a -x 
x 


dx + C =a In 


f 

a + 


V 


sfa 2 



sla 2 -x 2 +C 


und keine anderen dieser Differentialgleichung geniigen (CeR beliebig; siehe Formel Nr. 9 im An- 
hang 1). Wegen der Anfangsbedingung y(a) = 0 muss C = 0 sein. Die gesuchte Gleichung der Traktrix 
ist also 


y{x) - a In 


a + 


sla 2 


\~2 2 

sla -x . 


Was aber weder Claudius Perraltus noch Gothofredus Guilielmus Leibnitius im entfemtesten ahnen konnten: Wenn man ihre 
verspielte Traktrix um die j-Achse rotieren lasst, entsteht eine Flache mit konstanter negativer GauBscher Kriimmung - und 
diese „Pseudosphare“ ist nichts weniger als ein Modell fur die Nichteuklidische Geometrie Lob ats chew skis (1792-1856). 
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1. Attempt: with the 77-92 

Trying to find the tractrix by integration using the 77-92 you will find another result as given in 
Heuser’s book: 



' F** Y FMT Y FS Y Ffi^ Y ) 
Calc Other |Prgn I CMC lean Up| 

.11 -!•*-**. 


\{ X J 

Q. X 

a ln(lxl) - a* lnQa 2 - x 2 - aj -Ja 2 - x 2 

.[( J- 3 * -* 2 

cix | x > 0 

Jl X . 

-[a • lnlJa 2 - x 2 - aj - a ln(x) + Ja 2 - x 2 J 

x A 2>/x.x>lx>0 

MAIN RAD AUTO FUNC 2/50 



We integrate as usual. 

Integration with condition x > 0. 

The tractrix shall start in point (2,0) 



The graph shows the - unexpected - result that we find 
the curve reflected with respect to they-axis. Why that? 

Something must have gone wrong with the sign ... 


2. Attempt: again with the 77-92 (PLUS Module or 77-89 and the Voyage 200 as well) 

The PLUS (and the 77-89) offers the opportunity to solve differential equations using the deSol ve 

command. Applying this command to integrate yields the same result as before: 


m 

y PW Y Y l 

fl 1 gebra |Ca 1 c {Other |Prgn 1 0 1 

[dean Upf 

■ brapn "^•inu^-x — - t - 11 

riKX) + - X 


Done 


■deSolve[y‘ = ~ x2 , x, y] | x > Q 


y = “2- ln(j4 - x 2 - 2 ) + 2 ln(x) - J 4 - x 2 + £1 

■ Define yl(x) =y = "2 ln(j4 - x 2 - 2 ) + 2- lr* 

Done 





MftIN 

RAD AUTO 

FUNC fi/30 


P^l2^jTr F |cefReG7aphlH F a 5 tYl0» SW 1 





MAIN RAD HUTU FUNC 



You’d better choose a numerical approach: 


(Don’t forget first to set MODE 6:DIFF EQUA- 
TIONS) 
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If you take care to have a defined radicand you will obtain the correct “drag curve”. 



Let’s go on exploring ... 

3. Attempt: once more with the 77s 

If we want to find the reason for this “bug” then we have to investigate the expression to be integrated 
very accurately. Just with the first summand we are lucky: the root is defined only for |x| > 0 and re- 
turns a value between 0 and 2. Subtracting 2 leads to the fact that the whole summand is undefined for 
any real number. It is more than astonishing that a curve is plotted yet! 



Let us now manipulate this expression (unfortunately we can do this only “manually” - or as we Aus- 
trians are saying “handisch = zu FuB (by foot)” - because from my point of view the 77 is too poor in 
manipulating logarithmic expressions) and then plot the respective expression. By collecting the first 
two summands we change the domain and y2 ( x ) is now in fact defined for [-2,0]. 



Dragging the factor -2 into the log-function the domain changes again: y 3 ( x ) and y4 ( x ) are defined 
for [-2,0] u ]0,2]. 
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Plotting function y4(x) we obtain the same curve, what is quite clear. We have to remove the root 
from the denominator to receive the requested and expected result because the log-expression is now 
defined for only ]0,2]. 



We can see that the manipulations connected with logarithm, root and square are influencing the do- 
main of the graph in a very high degree. The more accurate inspection has shown that the “bug” in fact 
is no bug. But should we not wish to have immediately the solution as presented in y5 ( x ) . As a user 
only I cannot estimate how overbold this wish really is ... 

Now let us see how DERIVE is treating this integral: 


4. Attempt: with DERIVE 4.10 

Here manipulation of the expres- 
sion goes one step further and so it 
is defined for ]-2,0]. Function and 
graph are here - which is a differ- 
ence to the first 77-solution - corre- 
sponding. 

Again we need some manipula- 
tions (like on the TI but with the 
difference that the decisive step - 
removing the root from the de- 
nominator - is performed by 
DERIVE) to find the requested 
result (see #6:)) 
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5. Attempt: with Maple V (R4 Student Version) 


In this connection it might be in- 
teresting to compare how other 
CAS are dealing with this problem. 
Maple returns a new representation 
of the integral. But you will fail 
plotting the function graph - mes- 
sage of the system: “empty plot”. 
Obviously here is nothing to plot: 
the Area function arctanh(x) is 
defined for |x|<l. This leads to the 
consequence that for no x e [0,2] 
the argument lies in the domain of 
arctanh(x). 


STUDENT > expr} \- 


exprl '=-*Ja 2 -x 2 j 


STUDENT > stmplify(subs(d = 2, expr I)) 


* arctanh — — h 1 ■ “ t 
^Aa 2 J a 2 -x* 


4-x + 2 arctj 


[STUDENT > £1 := x->- (4-x~2) * (1/2) +2*arctanh (2/ (4-x A 2 ) * (1/2 ) ) ; 


fl :*=jc -> \J 4 -x + 2 arctanh 


STUDENT > plot(fl(x) ,x=0 . .2) ; 

Plotting *r«roir» *mpty plot 


1 Y -|- jc 

Using the relation arctanh(x) = — In (with | x |< 1) we find one of the known representations of our 

2 1 — x 

problem: 


V(4 - x ) 


V(4 - x ) = LN - 


2 2 
(V( 4 - x ) + 2) 


- V(4 - x ) 


V( 4 - X ) 


The log expression can again be transformed such that it represents the expected curve. 


6. Attempt: with Mathematica 3.0 (Teacher’s Edition) 



MATHEMATICA behaves very brilliant. 
Version 2.21 returns exactly the function 
given by Heusler, version 3.0 a slightly 
modified one. 

It leads directly to the graph. 


[1] H. Heuser, Gewohnliche Differentialgleichungen, Teubner Verlag, Stuttgart 1995 

[2] Thomas Weth, A Lexicon of Curves (5), DERIVE Newsletter #16 
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DIOPHANTINE EQUATIONS (2) 


Alfonso J. Poblacion Saez 
E.U.Politecnica, Valladolid, Spain 


Third Problem (presented in ACDC 6/DNL#29): The tiles problem could be solved as follows: 
let x 2 be the number of tiles used for the first nave, and y 2 for the third one. Then 


Putting y/2 = z, we have x 2 - 13z 2 = 1, which is Pell’s equation. 

For solving it we will use the second file explained in ACDC 7, adding some new functions: 

#10: [M (n) :=,L(n) :=] 

#11: M_AUX(n) :=IF(n=l, 0,w SUB 1) 

#12: L_AUX(n) : =IF (n=l , 1 , num A 2-w SUB 1 A 2) 

#13: L (n) :=IF(n>2, (num A 2-M(n) A 2) /L(n-l) ,L_AUX(n) ) 

#14: M(n) :=IF(n>2,L(n-l)*w SUB (n-1 ) -M (n-1 ) , M_AUX (n) ) 

#15: TABL (k) :=VECTOR ( [n, w SUB n, P (n) , Q (n) , M (n) , L (n) ] , n, 1 , k) ' 

Pell’s equation verifies that P(ft) 2 - num 2 Q(n) 2 = (-1)” L(«+l). In our case num = V 13, and you will 
obtain the data for finding the solution of our equation, adding for example: 

#16: S_L(k) : =VECTOR (IF ( P (n) A 2-num A 2 *Q (n) A 2=l, [P (n) , Q (n) ] , 0) , n, 1, k) 

#17 : SOL (k) : =SELECT (i/=0 f i, S_L (k) ) 

and then approximating SOL ( 1 0 ) . The final sum of tiles then is: 2x 2 = 2 * 649 2 = 842 402. 

For the last two problems we will make use of a trick involving the function ROW_REDUCE ( * ] . 

Fourth Problem: 




Then simplify 


10 2 c 


25 16 74 c 


R0W_R EDUCE 


52 -3S 47 c 


0 1 


2 


so that we can conclude that the two possibilities are 234 and 468. 


SOLUTION5([25 a + 16 b = 74 c, 52 a - BS b = 47 c], [a, b]) 


2 c 


does the same job, 
Josef. 


DERIVE' s' ROW_REDUCE function is equivalent to the 77’ s rref() function! 
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But the fifth problem is not so immediate to be solved. You have to write nine equations: eight of 
them for rows, columns and diagonals are equal to 15; the ninth with the sum of digits from 1 through 
9, which is 45. Then we can use ROW_REDUCE like in the former question. 

As the matrix is singular, some rows result in zeros. Hence we can write a, = 1 0 /. x 2 = 10-z, 

X 3 = -5+t+z, X 4 = - 1 0+2 t+z, a ' 5 = 5, Xf, = 20-2t, x 7 = 1 5 —t—z, x$ = z, and xg = t. Now we can calculate all 
possibilities which can created in a statement like 

VECTOR (VECTOR ( [10-t, 10-z , -5+t+z , -10+2 t+z , 20-2 t-z , 15-t-z , t, z ] ,z,l,9),t,l,9) 

and then select from it the solutions which have no numbers repeated (how?). 

(MAGSQ.MTH) 

Finally a challenge for which I have been thinking about a lot of time: Has anyone any idea 
how to obtain the 880 different order 4 magic squares with DERIVE ? 


I gave a workshop at TIME 2014 - “Brain Twisters - and how to tackle them”. There were 
also problems like Alfonso’s fourth problem, Josef. 

vars(v, k) creates all permutations (with repetitions) of order k of the elements given in vec- 
tor v. The two conditions are described as equations. SELECT finds out all three element vec- 
tors (out of 1000) which satisfy both conditions. 

vars(v, k, b, k_, m_ := 0, n_, s_ := [], t_) := 

Loop 

b := DIM(v) 

If m_ = b A k 

RETURN REVERSE (s_) 
k_ := k 
n_ := m_ 
t_ := [] 

Loop 

t_ := ADJOINfvj. (M0D(n_, b) +1), t_) 
n_ := FLOOR (n_, b) 
k_ :- 1 

If k_ = 0 exit 
s_ := ADJ0IN(t_, s_) 
m_ :+ 1 


all := vars([0, 1, 2, 3, 4, 5, 6, 7, 8, 9], 3) 

["a := v , b := v , c := v 
L 1 2 3. 


condl := 100 -a + 10-b + c = 26 -(a + b + c) 
cond2 := 100 -c + 10-b + a = 48- (a + b + c) 


SELECT(condl a cond2, v, all) 


" 0 0 0 ‘ 

2 3 4 

. 4 6 8 . 


needs 0.3 sec calculation time. 
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What you always wanted to know about the Undefined Variable Error 

but never dared to ask 

Wolfgang Propper Niimberg, Germany 

Each 11-92 user is familiar with the error message shown in the title. From the experience the under- 
signed has gained in many teacher training courses on the TI-92 (as well as from his own inadequacy) 
it might even be one of the most frequently occurring errors. In this contribution the reasons for it are 
identified and methods are shown, how one can avoid and overcome it. 

1 . As the error message says, the TI-92 has come across an undefined variable. Depending on the 
situation it can be easy or more difficult to identify this (or those) variablc(s). 

A first remark: This error message only occurs with graphing and/or tabulating functions (which 
are closely related). Generally algebraic manipulations (like limits, derivatives, etc.) in the Home 
Screen do not ask whether a variable is declared or not. 

When searching for the reason, one should look at the error message exactly: There are two 
slightly different forms of the message. 



With the message shown in Fig. 1, i.e. the one with two alternatives, you have already nearly 
found the reason: 

With an I ENTER 1 you return to the Home Screen and the entry line shows for example (s. Fig. 3) 



Looking into the Y= Editor (s. Fig. 4) belonging to it, one sees that the function which is said to be 
defined in the entry line is activated there. Possibly in this example the exponential function was 
to be plotted. And indeed it is called exp(x) in nearly every programming language. With the TI- 
92 however it is reached by I2ndl [LNl . 

The same case also occurs when one refers to a function in the Y= Editor, which was declared 
once, but has since been deleted (e.g. with DelVar f or by pressing (£6]). However the correspond- 
ing line in the Y= Editor is still activated. Changing to the Graph Screen causes the error message. 
The case described here in Graph Mode "Function'’ is analogous to the four other plot modes. 

A similar situation to the case described above occurs when declaring a term in x with an addi- 
tional parameter in theW= Editor (s. Fig, 5). When assigning values to that parameter using the 
"with" operator ([2nd] [k])> everything is okay and the TI-92 plots (like in this case) the graph of a 
family of curves (s. Fig. 6). 



However one must be careful: When declaring a function with parameters in an analogous way in 
the Home Screen (Fig. 7) and then trying to plot the family (Fig. 8), one gets "shipwrecked". The 
TT-97 does not pass the "with" statement to the symbolic term f(x). Instead it generates an error mes- 
sage like in Fig. 1 and displays it on I ENTER I like in Fig. 3. 
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■ Defin e f(x) =x^ 

rw»Mvni|P 


Fig. 7 


(Ti7^YFrrYT5 YfhYfs^Y re- Y r- r 

|»^|Zoom|Editp |fll 1 Istylejp* , 

.•.PLOTS 

^yl*f(x) 1 1 * C "4 0 5> 

a? 


Fig. 8 


(This flaw has been eliminated in the TI-92 Plus and the situation shown in Fig. 7/8 works with- 
out any error message.) 

There are two ways to overcome the errors described above: 

• One either deletes or changes the functions that caused the error message - either in the Home 
Screen or in the Y= Editor - or 

• one deactivates the functions with (£4] (in the Y= Editor) or executes the FnOff statement in the 
Home Screen fl£4] [0). 

A totally different way to generate the undefined variable error involves using the Graph statement 
of the TI-92. It is copied from the Other menu in the Home Screen to the entry line by fF4| f2l. 

Every body will regard it a truism that the statement Graph f(x) will cause an error message if the 
function f was not declared. The message comes in the Graph Screen and can only be closed by 
1ESC1 (Fig. 2). Subsequently one returns to the Home Screen with [T] (jL) an d finds, as a result of 
the Graph statement, the familiar error message (Fig. 9, 1 st line). 

The case becomes interesting if one now declares a function with another name and tries to plot it 
with Graph g(x). The same error message occurs. One feels totally lost if a direct Graph statement 
as shown in Fig. 9 is added. Here the same undefined variable error arises! 


■ Graph fO) Error: Undefined variable 

■Define g<x)=x 2 Done 
■Graph g(x) Error: Undefined variable 
■Graph x-1 Errors Undefined variable 


Undefined variable 


fESC=CRNCELj 


Fig. 9 


expl<x)=f <x> 

MAIM RAD AUTO 


Fig. 10 


An explanation is found, in the Ta ble Sc reen ( T»l (~y~|). Here an undefined variable errror message 
is issued first. After closing it with IESCI one sees the reason (Fig. 10). The table for an expression 
exp1(x)=f(x) with an undeclared function f is to be generated. That means the Graph statement has 
done more than its name says. It did not only try to plot in the Graph Screen but has also written 
its argument into the Table Screen. 

b, fV :^Tir=r- r-^rnl 3 ) 


I?WU& 

lW*3** W !■ 

'•wllm** 


X 
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G. 

Hl 
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l. 

0. 




2. 

4. 

1. 




3. 

m 

2. 




4. 


3. 




5. 

iL 

4. 





exp3<x>=l»l 

MAIN RAD AUTO 


Fig. 1 1 


Fig. 12 


After deleting the highlighted f(x) with the ICLEARI key (and not forgetting IEMTERI !) one finds two 
tables which are numbered 2 and 3 (Fig. 1 1). One can be easily convinced that these are the tables 
of the two functions which were to be plotted with the 2 nd and 3 rd Graph statement. In the Graph 
Screen (|TJ [r]) one sees that the plot statements have been executed (Fig. 12). 

The easiest way to remove the error message in this case is the ClrGraph statement (it can be 
found in [F4) GO). It deletes all entries in the Table Screen including the cause for the error. Ex- 
tremely cautious people use a ClrGraph prior to each Graph statement. 

3. Finally, a very insidious situation: 

The TI-92 giudebook deals on page 188 ff. with a "Preview of Statistics and Data Plots". Carry' 
out steps 1 to 4 and 12 to 14. Now switch with [♦] [r] to the Graph Screen. Normally you will 
see nothing because your coordinate system will not be dimensioned in the right way. A simple 
ZoomData (fF2l [9]) however should produce a picture like Fig. 13. One can be easily convinced, 
(for instance by tracing: (F3)) that the set of points is nothing but the data (pairs of numbers) stored 
in build. 
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If you now change to the Home Screen, delete the data variable with DelVar build, and then 
switch to the Graph Screen again with [*] l~R~l , our "popular" error message is back again. 
(Instead of deleting the variable build and possibly having to repeat the entry procedure mentioned 
above, one can rename it in the Var-Link Screen (|2ndl Q) using [FT) {3] there.) 


Of course one knows the reason for the error. However if this data variable was deleted some time 
previously and you have not accessed the Graph Screen since, you may not remember the dele- 
tion and have difficulty localizing the error. You will not find any illegal Graph statement in the 
Home Screen and you will find practically no hints to that error In the Y= Editor. A close exami- 
nation shows the word PLOTS with an up arrow (t) between the Y1= line and the toolbar in very 
small characters. A few © yield the picture shown in Fig. 14. The little tick in front of PLOT1 
indicates that this line is active. With the IH) key this line can be deactivated or the PLOT can be 
deletetd with ICLEAR1 . 


Conclusion: If the undefined variable error occurs after the ClrGraph statement has been executed 
and all functions in the Y= Editor have been deactivated, then there onla remains to examine the 
PLOTS area (which is the same in all graph modes). 


In case, however, you find other occurences of this error please let me know. My eMail-adress is: 
w-proepper@wpro.firanken.de 


The next contribution has its source in a lecture given by Dr. Maria Koth at the Vienna University at 
the occasion of a so called “Teachers’ Day”. Every year one day in spring time - the first Friday after 
Easter - the university opens its doors for math teachers. Professors and teachers give lectures on 
various fields of interests for secondary school teachers. Dr. Koth is one of the university teachers who 
try to inform the students about CAS and how to use it in school. Her lectures are always well pre- 
pared and I have the pleasure to present the essentials of her talk “Graphics with DERIVE”. I believe 
that it might be useful for beginners and for more advanced users of DERIVE as well. Josef 


Computergrafik mit DERIVE - Computer Graphics with DERIVE 

Dr. Maria Koth, University of Vienna, Austria 


Contents: 

0. Introduction 

1 . Elementary graphics with lines 

2 . Regular polygons 

3 . The Koch-Curve 

4. The Sierpinski-Triangle 

5. The Chaos Game 

0. Introduction 

First take care for the following settings: Work in Approximate Mode with not more than 3 digits to 
display. In the Plot Window switch off the Axes and Color Cycling. Point Size should be Small and 
Points should be Connected. 
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1. Elementary Graphics with Lines 

1 . 1 Graphics over a right angle 



We choose (0,0), (0,^) and («, 0) as endpoints 
of the two segments because then it is easy to 
shift, stretch or shrink the graphic using Cen- 
ter and Scale. 


The two segments on both axes with equal lengths are 
divided in n equal intervals. 

Every i th dividing point on one segment is connected 
with the (n-i) th dividing point on the other segment. 
Then the following graphic is created: 



To produce a similar graphic using a programming language needs a lot of work. The 64 DERIVE pro- 
gram” consists of only one single line: 

GRAFIKl (n) := VECTOR ( [i, 0; 0, n-i] , i, 0, n) 


Enter now GRAFIKl (20) and plot, done. 

(Later we will try to reproduce the graphics with Tl-Nspire , Josef) 

You can describe the graphics below by simple DERIVE commands very similar. 

GRAFIK2 (n) := VECTOR ( [ i , 0 ; 0 , n-i ; -i , 0 ; 0 , i-n; i , 0 ] , i,0,n) 




GRAFIK2 (20) SUB 16 


GRAFIK2 (20) 
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GRAFIK3 (n) := VECTOR ( [ 0 , n-i ; i, 0 ; n, i ] , i, 0, n) 




GRAFIK4 (n) 


VECTOR ( [0,n-i;i,0;n,i;n-i,n;0,n-i] , i , 0 , n) 




GRAFIK4 (30) 


1 .2 Graphics over a 60° angle 

Starting with two segments enclosing an angle of 60° the segments form two sides of an equilateral 
triangle. From lower secondary school the pupils (should) know the altitude of such a triangle of side 

1 t n^J 3 

length n with . 

GRAFIK5 (n) : = VECTOR ( [n-i , 0 ; 0 . 5 i,0.5 i V 3 ] A i , 0 , n) 
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The following patterns are produced by GRAF I K 6 (n) , GRAFIK7 (n) and GRAFIK8 (n) . 

GRAFIK6 (n) : = VECTOR ( [i, 0; 0 . 5 (n-i),0.5 (n-i) a/ 3; -0.5 i,0.5 i a/3; 

i-n,0;-0.5 i,-0.5 i V3;0.5 (n-i), 0.5 (i-n) V3; i, 0 ] , i, 0, n) 

GRAFIK7 (n) := VECTOR ( [0.5 (n-i), 0.5 (n-i) a/ 3 ; i, 0 ; n-0 . 5 i, 

0.5 i a/ 3] , i, 0, n) 

GRAFIK8 (n) : = VECTOR ( [ i , 0 ; n-0 . 5 i,0.5 i a/ 3;0.5 (n-i), 0.5 (n-i) a/3; 

-0.5 n+i,0.5 n V3;-0.5 i,0.5 i a/ 3;-0.5 (n+i) , 

0.5 (n-i) a/ 3; i-n, 0] , i, 0, n) 



GRAFIK6 (15) 


GRAFIK7 (40) 



Try this: 

GRAFIK9 (n) : = VECTOR ( [ i , 0 ; n-0 . 5 i,0.5 i a/3;0.5 (n-i), 0.5 (n-i) a/3; 

-0.5 n+i, 0.5 n a/3;— 0.5 i,0.5 i a/3 ; — 0 . 5 (n+i), 

0.5 (n-i) a/3; i-n, 0; -0 . 5 (n+i), 0.5 (i-n) V3;-0.5 i, 
-0.5 i a/3; -0.5 n+i, -0.5 n a/3;0.5 (n-i), 0.5 (i-n) a/ 3; 
n-0. 5 i,-0.5 i V3;i,0],i,0,n) 


GRAFIK9 (50) 
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2. Regular Polygons 

2 . 1 Plotting regular polygons 

The coordinates of a regular polygon with n vertices, its center in (0,0) and r as radius of its circumcir- 
cle are given by: 


( 2ni\ 

. ( 2 ni \ 

r-cos 

,r-sin 

n ) 

« )_ 


Using the DERIVE syntax this n-gon is described by the following list of points: 

POLYGON (n, r) == VECTOR ( [r -COS (27ii/n) , r • SIN (27ii/n) ] , i, 0, n) . 

There is another possibility to fix the vertices by: 

[r -cos(z),r -sin(/)] with 0 < / < 2/r step — . 

POLYGON2 (n, r) == VECTOR ( [r -COS (i) , r • SIN (i) ] , i, 0 , 271, 27i/n) . 

Try POLYGON (9,3) or POLYGON2 (12,4). Explore the VECTOR command with 
VECTOR (POLYGON2 (7,t),t,0.5,5,0.5) and VECTOR ( POLYGON2 ( t , 4 ) , t , 1 , 12 ) . 



The line graphs with TI-NspireCAS & LUA ( by Josef Bohm) 

Now - several years later - times of DERIVE, TI-92 and even Voyage 200 seem to have 
passed. Instead of this we have new tools like TI-NspireCAS, GeoGebra, CASIO, ... 

I will show how to create Maria Koth’s line graphs with TI-NspireCAS, which is not so easy 
because we miss my favourite DERIVE command - VECTOR. The students - and we too - are 
using vector as a tool, nothing else than a black box. So let’s provide another black box for 
Nspire. The important fact is finding out the “rule” for generating the petty families of seg- 
ments given by the matrices as first argument in the VECTOR construct. 

Another problem is that we cannot program the graph window directly - as it was possible 
with TI-92 and Voyage 200. So, how to plot the family of segments? 
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One solution is working with a LUA-script, and I will demonstrate how to do this: 

First of all I need a program 
for calculating the coordi- 
nates of the vertices of the 
segments in the right order. 

Then I open a Notes page 
and define the rules as ma- 
trices from grl through gr9 - 
and can find my own rules, 
too. 

For using the variables in 
the LUA-script it is neces- 
sary to work in approximate 
mode (i.e. enter at least one 
number in the “rule” as a 
decimal number). 


Igraph 

Define lgraph(rc//e,/j_)= 

Prgm 

Locator, k 

Ml]} 

£J:=dim{™fe)[l] 

For 7 , 0 f n_ 

For k f 1 ,d 

p/:=rafe[fr]|i=7_ and n\=n_ 
l. =au gment(au gmentf/, {pt[ 1, l] }), {p([ 1, 2 ] }) 
End For 
End For 
m\ =dirn(/) 

End Prgm 


6/10 


B 


grl: 


gr2:= 


gr3: 


gr4: 


gr5: 


7 

0 

► 

2 

0 


0 

72-2 


0 

72-2 


7 

0 


2 

0 


0 

72-2 


0 

l 72-2 


"7 

0 

► 

-2 0 


0 

2-72 


0 

l 2-72 


7 

0 


2 

0 


0 

72-2 


0 

72-2 


7 

0 


2 

0 


72 

2 


72 

2 


0 

72 - 

-2 


0 

72-2 

2 

0 


2 

0 




► 



72 

2 



72 

2 

72 - 

2 72 


72-2 

72 

0 

72 - 

-2 


0 

72-2 


- 


n-i 0 
0.5-7 0.5- v 

i 0 

w-0.5-1 0.5- v 

0.5- {n-i) 0.5- {n-i) m ^3 


n-i 0 

0.5-7 0.866025-7 


77 — 0 . 5*2 0.866025*7 

-0.5- (i-n) 0.866025 ■ (i-n) 




Then I insert a LUA-script. 

(For informing about LUA you may go to 

http://compasstech.com.au/TNS Authorinq/Scriptinq/index.html 


and explore Steve Arnold’s great LUA website.) 
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File Edit View 





hi Problem 1 



Insert 
Problem 
S Page 


Tools Window Help 


di Calculator 
Graphs 
Geometry 

'J Lists & Spreadsheet 
J3 Data & Statistics 
J Notes 

Vernier DataQuest™ 

3 Image... 

J Video... 

PublishView™ Text Box 

$ Program Editor 




D Insert Script 

r 

■ Edit Script [) 

1 r f. 


More LUA resources are given at the end of 
this article. 

Another LUA script is presented in the next 
contribution. 

(Steve recommends working with later ver- 
sion apilevel ‘2.2’.) 
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Run the program from the Notes page. 



o. j \tt i if cm? yrr rj 


U, J \l i UJ 

v* j \i ttj 

gr9:= 

i-n 0 

► 

i-n 

0 


0.5- (n+i) 0.5- (i-n)- 


-0.5- (i+n) 

0.866025- (i-n) 


0.5" i “0 .5- i m J 3 


-0.5 -i 

-0.866025- i 


-0.5 'n+i -0.5- n- J 3 


i-0.5- n 

-0.866025- n 


0.5 - (n-i) 0.5- (i-n)- Ji 


"0,5* (z'-h) 

0.866025- {i-n) 


n-0.5-i -0.5- r 


h-0.5 b i 

-0.866025- i 



i 

0 


i 0 




lgraph(gr9,5o) ► Done 





Just overwrite the parameters in the 
program call. Don’t enter it a second time. 

On the handheld we have to take n= 20. 
Other wise the graph would become too 
large. 


http://www.lua.org/ 
http://lua.gts-stolberg.de/ 
http ://www. t3vlaanderen.be/ fileadmin/ 13 -be/ cahiers/ cahier_3 5 .pdf 
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The Triangle of PASCAL 

Josef Bohm, Wiirmla, Austria 

On the last page I’d like to present a nice program which brought a lot of fun in some of my math les- 
sons. My students explored the binomial theorem, built up the Triangle of Pascal and learned how to 
use its numbers finding the expansion of powers of binomials considering additionally the pattern in 
its exponents. 

As a home work I gave them a sheet of paper with two 2 1 lines Pascal Triangles and asked them to 
mark in one figure all odd numbers and in the other one all numbers which are divisible by 3. 

They were very much surprised about the nice pattern which appeared under their hands. Then I 
transmitted a 77-92 program pas ( ) to their 77s and asked them first to remove the axes from the 
graph screen. As they had no plotting experience on the 771 need not caring for the [Y=] -editor. 

They ran pas(50,2,0,1 ) and then they started investigating: 


past n , a , r , fit) 

Prgm 

Local 

Define pt ( x , y )=Prgm 
PxlOn x,y:Pxl0n x,y:Px!0n x+l,y 
PxIOn x,y-l:PxlOn x+l,y-l 
EndPrgm 
Cl rDraw 
FnOff 
For i ,0,n 
For j.O^i 
nCr ( i , j )-s-z 

If mod ( z , a ) = r and et=l or 
mod(z,a)^r and 9t=0 Then 
pt( 1+2*1 , 117 -2*i+4*j ) 

Endlf 

EndFor 

EndFor 

FnOn 

EndPrgm 
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pas ( n , a , r , 0t ) gives the first n rows of the Triangle. The numbers which leave the remainder r by 
division by a are marked as a black spot (if 6 1 = 1). If 8 1 = 0 then all positions are marked where the 
remainders is not r (i.e. the negative picture). 

The students had the task to find nice PASCAL pictures in order to produce their PASCAL GAL- 
LERY. Another challenge was to find out the parameters of pas() for given pictures. (The third one 
was a bit tricky for them.) 

Gunter Scheu wrote an article on 46 Explorations in the Triangle of Pascal ” in DNL#9. Then - this was 
in 1993 - I was inspired by his contribution and by a chapter in Peitgen’s very recommendable book 
44 Fractals for the Classroom ” and provided the students with some DERIVE functions to produce simi- 
lar figures (see t ri pas . dfw) 

The latest versions of DERIVE offered the possibility for programming (see pascal _new. dfw). 



So we can do it with DERIVE and on the TI-92 but how to do it with TI-NspireCAS? 

We can again work with a LUA script, but we can also plot the respective points as a scatter 
diagram on a Geometry page or as a diagram in the statistics page. See the following screen 
shots. Base of all is a program similar to pas ( ) given above. 


pas 

Define pas {nn , aa , rr,s ) = 

Prgrn 

Local ij,z,txt 

MOWO} 

n\=nn\a\=aa\t:=tr 

ixt = { "is not remainder is remainder" } 

For i,Q,nn 
For y,0,/ 

z\ =mo(|n Cr(/ 5 /) , aa ) 

If z=tr and s=0 or z±tr and s=l Then 
It:=augmen{fr, {-4- /+8-j}): fy:=augment(fy,{ 6- /}) 
Endlf 
EndFor 
EndFor 
nh=dim(l\) 
fr:=fr/[s+l] 

EndPrgm 


We can run pas ( ) from the Notes: 


pas(90,3,0,l) * Done 

Divisor: 3 
0 is remainder 
Number of Points: 1372 
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You can insert the lists of coordinates in a Lists & Spread- 
sheet page. Please notice that I reversed the signs of the 
y-values because I wrote the program for being used in 
connection with the LUA-program. 


Running pas() with other arguments will change the plots 
immediately. 




A X_l 

ByJ 

- 

=lx 

=-ly 

1 

0 

0 

2 

-4 

-6 

3 

4 

-6 

4 

-8 

-12 

5 

0 

-12 

6 

8 

-12 

7 

-12 

-18 

8 

12 

-18 

9 

-16 

-24 

10 

-8 

-24 

11 

8 

-24 

19 

i a 

-9 A 


But you can also use the two lists 
to be presented in a Data & Sta- 
tistics page: pas(90,3,0,0) 

See finally the LUA-generated 
plot in a Geometry page. 


Scatter plot and Data plot as well 
are restricted to 2500 points. 

The LUA program allows more 
points to be plotted. 
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Many Thanks and Credits: 

I am very indebted to Steve Arnold who informed about a couple of improvements in my LUA- 
program due to the recent apilevel 2.2. Once more I am recommending visiting Steve’s website. You 
will find numerous Nspire applications and a very extended tutorial for LUA programs. 


SB, Script Editor 

File Edit Debug 

View 

Help 






® 


& 

m. 





a 

Suspend 

Resume 

Set 

Focus 

Step 

Step 

Enable 

Disable 

Set 

Script 

Script 

Script 

Script 

Into 

Over 

Breakpoints 

Breakpoints 

Permissions 


tripas 


1 

2 platform. apilevel = '2.2' 

3 local nl = 0; local lx = {};local ly = {}; function on. construction ( ) 

4 h=platf orm. window: height ( ) 

5 w=platf orm. window: width ( ) 

6 var . monitor ( rr nl rr ) ; var . monitor ( rr lx rr ) ; var . monitor ( rr ly rr ) 

7 end 

3 

9 — Import the variables 

10 function on. varChange ( ) 

11 nl= (var. recall ( rr nl rr ) ) 

12 lx=( var. recall ( rr lx rr ) ) ; ly= (var . recall ( rr ly rr ) ) 

13 platform. window: invalidate ( ) 

14 end 

15 

16 function on. resize (width, height) 

17 h=height; w= width 

13 nl= (var . recall ( rr nl rr ) ) 

19 — Impart the lists of coordinates 

20 lx= (var . recall ( rr lx rr ) ) ; ly= (var . recall ( rr ly rr ) ) 

21 end 

22 

23 function on. enterKey ( ) 

24platf orm. window: invalidate ( ) 

25 end 

26 lx = { } ; ly = { } 

27 for i=l,nl do 

28 lx[i] = 0 ; ly[i] = 0 

29 end 

30 function on. paint (gc) 

31 — Hotting- the Points 

32 for i=l,nl do 

33 gc: fillArc [w/2+lx[i] , 10+ly[i] ,5,5,0,360) 

34 end 

35 end 


http://compasstech.com.au/TNS Authorinq/Scriptinq/index.html 


p 64 


Josef Bohm: Patterns in PASCAL’S Triangle 


D-N-L#32 


Plots created with pascal _new. dfw: 
pasptvl(100, 2, 0) 


pasptvlCLSO, 6, 2) 



pasptv2(150 , 9, 0) 


W 

ww 



pasinvptv2(150 , 9, 0) 



